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We report a measurement of high-pr inclusive and direct photon production in p + p and 

d + Au collisions at ^/snn = 200 GeV at midrapidity (0 < 77 < 1). Photons from the decay 7r° — > 77 
were detected in the Barrel Electromagnetic Calorimeter of the STAR experiment at the Relativistic 
Heavy Ion Collider. The rj — >^ 77 decay was also observed and constituted the first r? measurement 
by STAR. The first direct photon cross section measurement by STAR is also presented, the 
signal was extracted statistically by subtracting the tt", rj, and a;(782) decay background from the 
inclusive photon distribution observed in the calorimeter. The analysis is described in detail, and 
the results are found to be in good agreement with earlier measurements and with next-to-leading 
order perturbative QCD calculations. 



PACS numbers: 13.85.Ni, 13.85.Qk, 13.87.Fh, 25.75.-q, 25.75.Dw 
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I. INTRODUCTION 

The high center-of-mass energy {y/sNN = 200 GcV) of 
the Relativistic Heavy Ion Collider (RHIC) opens up the 
hard scattering regime, which is accessed by measuring 
particle production at high transverse momentum px- 
The high-pT particles {px 3 GeV/c) originate from the 
fragmentation of partons that have scattered in the early 
stage of the collisions. Hence, in heavy-ion collisions 
the high-pT particles can be used to probe the produced 
medium of strongly interacting matter. A significant 
suppression of high-pT hadron production relative to a 
simple binary collision scaling from p+p has been observed 
at RHIC in central An -I- An collisions IJ. Furthermore, 
it was found that jet-like correlations opposite to trigger 
jets are suppressed, and that the azimuthal anisotropy 
in hadron emission persists out to very high pr [2H3]. 
In contrast, no suppression effects were seen in d + Au 
collisions [SMSj, which has led to the conclusion that the 
observations made in Au -I- Au are due to the high-density 
medium produced in such collisions and not to initial 
state effects. The most probable explanation to date is 
that the suppression is due to parton energy loss from 
induced gluon radiation (jet quenching) in the extremely 
hot and dense medium [9 . To quantitatively understand 
this behavior and, in particular, to separate hot from 
cold nuclear matter effects, such as Cronin effect [101 
and parton shadowing and antishadowing |11H13| , precise 
measurements of identified hadrons at high px in p + p 
and d + Au collisions are required [H] . 

Prompt photons have long been proposed as a powerful 
tool for studying the jet quenching via photon-jet cor- 
relations |15j . In the dominant hard photon production 
processes (quark-gluon Compton scattering and quark- 
antiquark annihilation) , the outgoing photon balances the 
momentum of its partner parton and has large enough 
mean free path to escape the collision system, provid- 
ing a calibrated probe for studying the energy loss and 
mean free path of the parton in the medium. In addition, 
prompt photons constitute a background for measuring 
the medium-induced production of photons in response 
to the energy deposited by that parton [16]. 

The thermal photon spectrum is directly related to 
the temperature of the hot and dense medium created 
in the heavy-ion collision, provided that it is in thermal 
equilibrium |17j . The measurement of such a spectrum 
requires a knowledge of the prompt photon background, 
which can be measured in p + p and d + Au systems 
that share multiple sources of photons with heavy-ion 
collisions but are not expected to produce an extended 
thermal system. 

The measurements of tt'^'s and direct photons in p + p 
collisions are also of specific interest for studies of the 
proton spin structure (see, e.g., Ref. [18.), which are un- 
derway at RHIC. A main objective of the RHIC spin 
program is to constrain the polarization of the gluons 
inside the proton [AG'(a;)]. The unpolarized cross sections 
provide a test of the next-to-leading order perturbative 
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QCD (NLO pQCD) framework, which is used to inter- 
pret the measured spin-dependent observables. 

In this paper, we present the first results for the high-pT 
tt'^, 77, and direct photon production in p + p and d + Au 
colhsions at ^/sNN = 200 GeV in the pseudorapidity range 
< ?7 < 1, measured by the STAR experiment at RHIC 
(except the cross section for 7r° production in p + p col- 
hsions, first presented in Ref. [E]). The STAR Barrel 
Electromagnetic Calorimeter was used to detect high-pT 
tt" and 77 mesons via their 77 decays. The direct photon 
signal was extracted statistically by subtracting the 7r°, 77, 
and w(782) decay background from the inclusive photon 
distribution observed in the calorimeter. The presented 
data constitute a necessary baseline for the measurements 
of TT*^, 77, and direct photon production in heavy-ion colli- 
sions at RHIC. Inclusive tt^ production was previously 
measured in STAR for low pT at midrapidity in Au -I- Au 
collisions at = 130 GeV [ID] and 200 GeV [H], and 

at the forward rapidities in p + p and d + Au collisions 
at ^/sj^ = 200 GeV [32]. STAR has also measured the 
production of other identified particles, such as tt*, 
p/p, and hadronic resonances [Till^gH^ . The PHENIX 
experiment at RHIC has also measured the 7r°, ij, and 
direct photon production at ^/snn = 200 GeV in a variety 
of colhsion systems, including p -I- p and d-l-Au [61 I26H29] . 

The paper is organized as follows. In section |TT| we 
describe the detectors that were used in this analysis. In 



section HI, we describe the data processing chain used 
to reconstruct photon candidates in the raw data. Sec- 
tions |IV] and |V] show how these photon candidates were 
used to calculate the yields of 7r° and r/, and direct pho- 



tons, respectively. Finally, in section VI we present the 
results and compare our data to the theoretical calcula- 
tions and to the measurements by other experiments. 



II. EXPERIMENTAL SETUP 

The STAR detector (Solenoidal Tracker At RHIC) ^ 
was designed primarily for measurements of hadron pro- 
duction in heavy-ion and proton-proton collisions over 
a large solid angle. For this purpose, tracking detectors 
with large acceptance and high granularity were placed 
inside a large- volume solenoidal magnetic field (0.5 T). 
The detector subsystems relevant for the present analysis 
are briefly described in the following sections. 



A. Time Projection Chamber 

The Time Projection Chamber (TPC) [H] is the central 
tracking device in STAR. It allows one to track charged 
particles, measure their momenta, and identify the particle 
species by measuring the ionization energy loss dE/dx. 

The TPC barrel measures 4.2 m in length, and has an 
inner radius of 0.5 m and an outer radius of 2 m. The 
TPC acceptance covers ±1.8 units in pseudorapidity and 
full azimuth. Particle momentum is measured in the range 



0.1-30 GeV/c. In this analysis, TPC tracks were used 
to reconstruct the interaction vertex and to identify the 
energy deposits of charged particles in the calorimeter. 



B. Forward TPC modules 

Two Forward Time Projection Chambers (FTPC) [32] 
extend the STAR tracking capability to the pseudora- 
pidity range 2.5 < I77I < 4. Each FTPC is a cylindrical 
volume with a diameter of 75 cm and a length of 120 cm, 
with radial drift field and pad readout chambers mounted 
on the outer cylindrical surface. Two such detectors were 
installed partially inside the main TPC, on both sides 
of the interaction point. In this analysis, the forward 
charged-track multiplicity recorded in the FTPC in the 
gold beam direction served as a measure of the centrality 
in d -I- Au collisions. 



C. Barrel Electromagnetic Calorimeter 

A Barrel Electromagnetic Calorimeter (BEMC) [55] 
was incrementally added to the STAR setup in 2001-2005 
to measure the energy deposited by high-pT^ photons and 
electrons and to provide a trigger signal. The calorimeter 
is located inside the magnet coil and surrounds the TPC, 
covering a pseudorapidity range j??! < 1 and full azimuth. 

The full calorimeter consists of two contiguous half- 
barrels, located east and west of the nominal interaction 
point, each of which is azimuthally segmented into 60 
modules. Each module is approximately 26 cm wide and 
covers 6° (105 mrad) in azimuth and one unit in pseudo- 
rapidity. The active depth is 23.5 cm, to which 6.6 cm of 
structural elements are added at the outer radius. Results 
presented in this paper used only the west calorimeter 
half-barrel (0 < 77 < 1), which was fully installed and 
calibrated in 2003-2005. 

The modules are segmented into 40 projective towers 
of lead-scintillator stacks, 2 in 1^9 and 20 in 77 direction. A 
tower covers 0.05 rad in Aip and 0.05 units in A77. Each 
calorimeter half-barrel is thus segmented into a total of 
2400 towers. Each tower consists of a stack of 20 layers 
of lead and 21 layers of scintillator. All these layers are 
5 mm thick, except the first two scintillator layers, which 
are 6 mm thick. A separate readout of these two layers 
provides the calorimeter preshower signal, which was not 
used in this analysis. A Shower Maximum Detector (see 
below) is positioned behind the fifth scintillator layer. The 
whole stack is held together by mechanical compression 
and friction between layers. From layer-by-layer tests 
of the BEMC optical system, together with an analysis 
of cosmic ray and beam test data, the nominal energy 
resolution of the calorimeter is estimated to be 5E/E = 

14%//E(GeV)® 1-5% M- 
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D. Shower Maximum Detector 

The Shower Maximum Detector (SMD) is a multi-wire 
proportional counter with strip readout. It is located at 
a depth of approximately 5.6 radiation lengths at 77 = 0, 
increasing to 7.9 radiation lengths at 77 = 1, including 
all material immediately in front of the calorimeter. The 
purpose of the SMD is to improve the spatial resolution of 
the calorimeter and to measure the shower profile. This is 
necessary because the transverse dimension of each tower 
(« 10 X 10 cm^) is much larger than the lateral spread of 
an electromagnetic shower. The improved resolution is 
essential to separate the two photon showers originating 
from the decay of high-momentum tt" and rj mesons. 

Independent cathode planes with strips along rj and (p 
directions allow the reconstruction of two projections of 
a shower. The coverage is A77 x A(p — 0.0064 x 0.1 rad 
for the 77-strips and 0.1 x 0.0064 rad for the (^-strips, each 
group of 2 X 2 towers covers 15 strips in each SMD plane 
behind it. In total, SMD contains 36000 strips. 



E. Trigger detectors 

In addition to the STAR barrel detectors, sampling 
liadronic calorimeters were placed at a distance of 18 m 
from the interaction point, on both sides of the experi- 
mental hall. In heavy-ion collisions, these Zero Degree 
Calorimeters (ZDC) |311[3S] measure the total energy of 
the unbound neutrons emitted from the nuclear fragments 
after a collision. The charged fragments of the collision are 
bent away by the RHIC dipole magnets upstream of the 
ZDCs. For the d+Au data used in this analysis, the ZDC 
provided a collision trigger by requiring the detection of 
at least one neutron in the gold beam direction. 

To provide a collision trigger in p + p collisions, Beam- 
Beam Counters (BBC) [33 |3Zj were mounted around the 
beam pipe beyond both poletips of the STAR magnet 
at a distance of 3.7 m from the interaction point. The 
detector consists of two sets of small and large hexagonal 
scintillator tiles arranged into a ring that covers pseudo- 
rapidities between 2.1 and 5.0. The minimum bias trigger 
required a coincidence of signals in at least one of the 18 
small BBC tiles on each side of the interaction region. 

The two BBC counters record timing signals that can 
be used to determine the time of flight for the forward 
fragments. The difference between these two flight times 
provides a measurement of the z position of the interaction 
vertex (zvcrt) to an accuracy of about 40 cm |38| . Events 
with large values of the time-of-flight difference, which 
indicate the passage of beam background, were rejected 
at the trigger level. The BBCs also served to measure 
the beam luminosity in p + p runs. 



III. DATA RECONSTRUCTION 
A. Datasets and statistics 

The data used in this analysis were taken in the d + Au 
run of 2003 and in the p + p run of 2005, both at ^/sNN = 
200 GeV. The integrated luminosity was 0.66 pb"^ for the 
p + p data and the equivalent nucleon-nucleon luminosity 
was 0.22 pb~^ for the d + Au data. The following trigger 
conditions were used: 

Minimum bias (MinBias) trigger in d + Au collisions 
This condition required the presence of at least one neu- 
tron signal in the ZDC in the gold beam direction. As 
determined from detailed simulations of the ZDC accep- 
tance [5], this trigger captured (95 ± 3)% of the total 
d + Au hadronic cross section of cTh^ck" = 2.21 ± 0.09 b. 

MinBias trigger in p + p collisions 
This condition required the coincidence of signals from two 
BBC tiles on opposing sides of the interaction point. Due 
to the dual-arm configuration, this trigger was sensitive 
to the non-singly diffractive (NSD) cross section, which 
is a sum of the non-diffractive and doubly diffractive cross 
sections. The total inelastic cross section is a sum of the 
NSD and singly diffractive cross sections. A MinBias 
cross section of (Tbbc = 26.1 ± 0.2 (stat.) ± 1.8 (syst.) mb 
was independently measured via van der Meer scans in 
dedicated accelerator runs [31]. This trigger captured 
(87±8)% of the p+p NSD cross section, as was determined 
from a detailed simulation of the BBC acceptance [T]. 
Correcting the BBC cross section for the acceptance, we 
obtained the NSD cross section cr^g^p = 30.0 ± 3.5 mb. 

HighTower trigger 
This condition required a transverse energy deposit Et 
above a predefined threshold in at least one calorimeter 
tower, in addition to satisfying the MinBias condition. 
This trigger enriched the recorded dataset with events that 
had a large Et- Two different thresholds were applied, 
defining the HighTower-1 and IIighTower-2 datasets. The 
nominal values of these thresholds were set to 2.6 and 
3.5 GeV in p + p, and to 2.5 and 4.5 GeV in d -I- Au 
runs. Prior to each run, all towers were equalized to give 
the uniform transverse energy response, by adjusting the 
high-voltage settings of the individual photoniultipliers. 

HighTower software filter 
The HighTower-triggered data were additionally filtered 
using a software implementation of the HighTower trigger. 
In this filter, the highest tower ADC value found in 
the event was required to exceed the same HighTower-l 
(IIighTower-2) threshold as the one that was used during 
the run. This filter was needed to remove events that were 
falsely triggered due to the presence of noisy channels 
(hot towers). Such channels were identified offline in a 
separate analysis and recorded in a database. In addition, 
the highest calibrated transverse energy of a tower in the 
event was required to exceed slightly higher thresholds 
{Et + 0-5 GeV) than those used during the run, to account 
for possible inaccuracy of the online calibration of the 



6 



xlO-' 



^ 2.0 
H 1-5 
1.0 
0.5 



(a) d + Au HighTower-2 
Beam bg cut 




0.2 0.4 0.6 0.8 



' Ptpo) 




FIG. 1: (a) Distribution of r = i5BEMc/(-E'BEMC + Ptpc) in 
d + Au events, which shows beam background at r > 0.8. The 
curve corresponds to a second order polynomial fit, constrained 
to pass through zero at r = 1, used to estimate the false 
rejection rate, (b) Distribution of r in p + p events. 



cause of the BBC coincidence requirement in the trigger 
and the timing cut on the BBC vertex position. Therefore, 
no cut on r was appHed to the p + p data. 

The residual beam background contamination in the 
d + Au MinBias trigger was estimated from an analysis of 
the empty RHIC bunches to be (5±1)% [14]. To estimate 
the residual background in our data, we analyzed a sample 
of 3 X 10^ MinBias tri ggers from unpaired RHIC bunches. 
These events were passed through the same reconstruc- 
tion procedure as other data. We observed that «10% of 
the fake triggers passed all cuts, and that none of these 
contained a reconstructed tt". The residual beam back- 
ground contamination in the tt*' yield was thus estimated 
to be 0.1 X 5% = 0.5% and considered to be neghgible. 



towers. This software filter also served to make the trigger 
efficiency for the Monte Carlo detector simulation and for 
the real data as close as possible. 

B. Beam background rejection 

During the data taking in 2003-2005, interactions of 
beam ions with material approximately 40 m upstream 
from the interaction region gave rise to particles that 
traversed the detector almost parallel to the beam di- 
rection. This source of background was eliminated by 
installing additional shielding in the RHIC tunnel for the 
subsequent runs. 

To identify events containing such background tracks, 
the ratio 

.EbeMC /in 

r = — (1) 

-Ebemc + Ptpc 

was calculated, where E'bemc is the total transverse en- 
ergy recorded in the BEMC and Ptpc is the transverse 
momentum sum of all charged tracks reconstructed in the 
TPC. In events containing background, r was large (close 
to 1) because photons from these background events de- 
posited a large amount of energy in the calorimeter, while 
the accompanying charged tracks were not reconstructed 
in the TPC, because they did not point to the vertex. 
Figure [l] shows the distributions of r for the d + Au and 
p-\-p data. The peak near unity in panel (a) indicates the 
presence of beam background in c?-|- Au collisions. Events 
with r > 0.8 were removed from the d-l- Au analysis. This 
cut rejected 3.4% of MinBias and 13% of HighTower-2 
events. From a polynomial fit to the d -\- Au distribution 
in the region r = 0.6-0.8 [curve in Fig. [ija)], the false 
rejection rate was estimated to be 3.6% in the d + Au 
HighTower-2 data and less than 1% in the other datasets. 
By studying this rejection rate as a function of i?BEMC; 
we estimated the potential distortion of the tt°, ij, and 
photon spectra due to the removal of these events to be 
below 1% in all datasets. 

Figure [ijb) shows the distribution of r for the 
HighTower p + p data. The background was negligible be- 



C. Determination of centralities 

To measure the centrality in d + Au collisions, we used 
the correlation between the impact parameter of the col- 
lision and the charged-track multiplicity in the forward 
direction. This correlation was established from a Monte 
Carlo Glauber simulation ^40r,-42j using, as an input, the 
Woods-Saxon nuclear matter density for the gold ion [33] 
and the Hulthen wave function of the deuteron [33] . In 
this simulation, the inelastic cross section for a nucleon- 
nucleon collision was taken to be ct;^^ = 42 mb. The pro- 
duced particles were propagated through a full GEANT [3S] 
simulation of the STAR detector. Both the charged track 
multiplicity and the number of nucleon-nucleon collisions 
simulated by the event generator were recorded. 

For the event-by-event centrality determination, we 
measured the multiplicity of tracks reconstructed in the 
FTPC module in the gold beam direction (iVFTPc)- Cen- 
trality bins were defined following the scheme used in 
other STAR publications [5]. The following quality cuts 
were applied to the reconstructed tracks: (i) at least 6 hits 
were required on the track; (ii) < 3 GeV/c, which guar- 
anteed that the track was fully contained in the FTPC 
acceptance; and (iii) distance of closest approach to the 
vertex had to be less than 3 cm. The multiplicity distribu- 
tions obtained from the d + Au data are shown in Fig. [2] 
for the MinBias, HighTower-1, and HighTower-2 triggers. 

Based on iVpTPCi the events were separated into three 
centrality classes: 0-20% most central, 20-40% mid- 
central, and 40-100% most peripheral, as indicated by 
the vertical lines in Fig. [2] Table |T] lists the A^ftpc ranges 
and the corresponding mean numbers of binary collisions 
((A'^coii)) obtained from the Glauber model, for each cen- 
trality class. The systematic uncertainties on {Ncou) were 
estimated by varying the Glauber model parameters. 



D. Vertex finding efficiency 

In p + p data, a vertex was reconstructed based on 
the tracking information for 65% of the MinBias events. 
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FIG. 2: Centrality selection in the d + An data, based on 
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defined, containing 0-20% most central, 20-40% mid-central, 
and 40-100% most peripheral events, respectively. 
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FIG. 3: Electron energy measured in the BEMC after calibra- 
tion, divided by the momentum measured in the TPC, in the 
p + p data. Solid line is a Gaussian fit, which shows that the 
peak is centered at unity. 



consequence, the TPC tracking efficiency is reduced for 
vertices located far from the center of the detector. 



TABLE I: Centrality classes defined for the d+Au data and the 
corresponding (Ncou) values [5]. The errors given for (Ncou) 
indicate the systematic uncertainty. 



Centrality class 


A'ftpc range 


(iVcoll) 


d -1- Au MinBias 




7.5 ±0.4 


0-20% most central 


> 17 


15. Oil. 1 


20-40% mid-central 


10-16 


10.2±1.0 


40-100% most peripheral 


< 10 


4.0 ±0.3 


p + p 




1 



For the remaining events, the vertex position in z was 
determined using the time information from the BBCs. 

In the d + Au HighTower data, the charged track multi- 
plicities were large enough to always have a reconstructed 
vertex. However, a vertex was missing in about 7% of 
the MinBias events, and cannot be recovered from BBC 
information because the BBC was not included in the 
d + Au MinBias trigger. Events without a vertex have 
low charged track multiplicity, and the contribution from 
these events to the 7r° yield above 1 GeV was assumed 
to be negligible [1^. Therefore, a correction for vertex 
inefficiency was applied as a constant normalization factor 
to the yield and its uncertainty contributed to the total 
normalization uncertainty of the measured cross sections. 

The vertex reconstruction efficiency in triggered d + Au 
MinBias events was Evert = 0.93 ± 0.01 However, 
this efficiency depends on the collision centrality, and we 
assumed that it was 100% for central events. Scaling the 
efficiency above by the ratio of peripheral to total number 
oi d + Au events, we obtained an efficiency correction 
factor of 0.88 ± 0.02 for the sample of peripheral events. 

Events with jzvert | > 60 cm were rejected in the analysis, 
because the amount of material traversed by a particle 
increases dramatically at large values of Invert |- As a 



E. Energy calibration of the calorimeter 

In the first step of the calorimeter calibration, the gains 
of the individual towers were matched to achieve an over- 
all uniform response of the detector. For this purpose, 
minimum ionizing particles (MIP) were used, by selecting 
the TPC tracks of sufficiently large momentum (greater 
than wl GeV/c). These tracks were extrapolated to the 
BEMC and the response spectra were accumulated, pro- 
vided that the track extrapolation was contained within 
one tower and that the track was isolated. For ci -I- Au data, 
the isolation criterion meant that no other tracks were 
found in a 3 X 3 patch around the tower; for p + p data, 
these neighboring towers were required to have no signal 
above noise. The peak positions of such MIP signals 
were used to calculate the tower-by-tower gain corrections 
needed to equalize the detector response [17]. 

In the second step, the energy scale was determined by 
comparing the momenta p of identified electrons in the 
TPC with the energies E recorded in the BEMC using 
the relation E = p ioi ultra-relativistic electrons. Figure [3] 
shows the distribution of E/p for a selected sample of at 
least 90% pure electrons in the p+p data at p > 2 GeV/c, 
after the calibration has been performed. The Gaussian fit 
to the central part of the electron peak demonstrates that 
the mean has been placed at unity. From a variation of 
the peak position with p, the systematic uncertainty of the 
electron calibration was conservatively estimated to be 5%. 
Within the present statistics, that calibration covers the 
momentum range only up to p = 6 GeV/c. Because the 
peak position is close to unity atp> 3.5 GeV/c, we assume 
that the assigned systematic uncertainty covers possible 
non-linearities at higher photon momenta p ^ 15 GeV/c 
probed in the present measurements. 

This calibration method takes advantage of the well 
understood TPC detector for the precise measurement 
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TABLE II: Cluster finder threshold values used in the analysis. 
Esccd, Ea_dd, and -Emin are the minimal energies for the seed 
hits, regular hits, and entire clusters, respectively. A^max is the 
maximal cluster size. 



Detector 


£.ccd (GeV) 


Sadd (GeV) 


Smin (GeV) 


max 


Towers 


0.35 


0.035 


0.02 


4 


SMD 


0.2 


0.0005 


0.1 


5 



of tlie electron track momentum in a wide range. A 
disadvantage is that it takes large statistics to calibrate 
the high-energy part of the spectrum. For this reason, 
only one global calibration constant was obtained. It was 
found that the current calibration is less reliable at the 
edges of the half-barrel. Therefore, the signals from the 
two 7y-rings at each side were removed from the analysis. 

The absolute energy calibration of the SMD was de- 
termined using the beam test data to an accuracy of 
about 20%. This analysis is not very sensitive to the ab- 
solute energy scale of the SMD, because the main energy 
mesurement was done with the towers. 



F. Particle reconstruction in the BEMC 

The first step in the photon reconstruction was to find 
clusters of energy deposits in the calorimeter by grouping 
adjacent hits that were likely to have originated from a 
single incident particle. The cluster finding algorithm was 
applied to the signals from BEMC towers and from each 
of the two SMD layers. 

The clustering started from the most energetic hit (seed) 
in a module and added neighboring hits of decreasing en- 
ergy to the cluster, until either a predefined maximal 
cluster size or a lower hit energy threshold were reached. 
The algorithm then proceeded to process the next seed. 
The threshold values are listed in Table [11} By construc- 
tion, the clusters were confined within a module and could 
not be shared by adjacent modules. However, the like- 
lihood of shower sharing between modules is considered 
to be low, since the modules are physically separated by 
«12 mm wide air gaps. The rj-ip position of each clus- 
ter was calculated as the energy-weighted mean of the 
individual hit positions within the cluster. 

After the tower and SMD clusters were found, they were 
combined into BEMC points, which closely corresponded 
to the impact points and energy deposits of particles 
that traversed the calorimeter. The procedure for form- 
ing the BEMC points is described in detail in Ref. |38j . 
The SMD information was essential because the mini- 
mal opening angle of the decay photons decreases with 
increasing energy of the parent n'^ . The spatial resolution 
of the BEMC towers alone is not sufficient to efficiently 
resolve the decay photons of tt^'s with p > 5 GeV/c. For 
this reason, only the BEMC points that contained tower, 
SMD-?7, and SMD-(/3 clusters were kept for the further 




FIG. 4: Distribution of the distances D between BEMC points 
and their closest tracks, obtained from p + p HighTower-1 data 
in the bin 4 < pt < 5 GeV/c. The curve shows a fit to Eq. ([sjl, 
the vertical line indicates the CPV cut. 

analysis of the High Tower data. In the analysis of MinBias 
data, used to obtain the 7r° signal at pt < 4 GeV/c, all 
reconstructed BEMC points were used, even when they 
did not contain SMD clusters. 

The SMD efficiency decreases rapidly and its energy 
resolution becomes poor with decreasing energy of the 
traversing particle, leading to significant fluctuations in 
the strip readout for E 2 GeV. Therefore, in the 
HighTower-1 data the SMD clusters were accepted only 
when they contained signals from at least two strips. This 
cut rejected a large fraction of the distorted and falsely 
spht SMD clusters, and reduced a possible effect of poor 
SMD response simulation at low energies. 

G. Charged particle veto using TPC 

A charged particle veto (CPV) cut was applied to reject 
the charged hadrons that were detected in the calorimeter. 
A charged hadron was recognized as a BEMC cluster with 
a TPC track pointing to it. The cluster was rejected if 
the distance D between the BEMC point and the closest 
TPC track in the r]-tp coordinates was 

D = + (A(^)2 < 0.04. (2) 

When a track was projected to the calorimeter surface, 
at a radius R = 220 cm, this cut corresponded to a linear 
separation RD « 10 cm in the pseudorapidity range of 
this measurement. The efficiency of this cut was 35% 
in the MinBias data and 71% in the HighTower-1 aird 
HighTower-2 data, in all p -I- p and d + Au datasets. The 
BEMC points remaining after this cut were considered 
to be photon candidates, and were combined into pairs, 
defining the set of tt" candidates. 

This veto introduced a false rejection of photon clus- 
ters if an unrelated charged particle happened to hit the 
calorimeter close to the cluster. Figure [4] shows the distri- 
bution of D observed in the p + p data. In this plot, one 
distinguishes the peak of real charged particles at small 
distances, superimposed on a random component seen 
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as a shoulder at larger distances. Assuming a uniform 
distribution of track projections in 77 and ip around the 
BEMC point, the radial distribution is given by 



f{D)^De-°''. 



(3) 



Here p is the charged track density in the vicinity of the 
photon. This parameter was obtained from a simultane- 
ous fit to the data in all bins of the event multiplicity M 
measured in the TPC, assuming a linear dependence on 
M, p — a + bM. The parametrization given by Eq. (|3| de- 
scribes well the random component, as shown by the curve 
in Fig. [4] The relative number of random coincidences 
that were falsely rejected was obtained by integrating the 
fitted curve up to the distance cut and weighting with the 
multiplicity distribution observed in each px bin. The 
resulting correction factor was e^p^ = 0.94 ± 0.02 for the 
p + p data and 0.89 ± 0.02 for the d -I- Au data. 

In the direct photon analysis, the purity of the photon 
candidate sample was more important than in the tt" 
analysis, therefore, a stronger cut RD < 15 cm was used. 
The correction factors were calculated to be e2pv — 0.95 ± 
0.02 for p + p and 0.93 ± 0.02 for d + Au data. The 
residual contamination by charged particles (C±) was 
estimated from the integrated excess of the D distribution 
over the fit to the random associations in the interval 
15 < RD < 25 cm, and was less than 5% for all pT bins. 

The uncertainties of these corrections contributed to a 
PT-independent systematic uncertainty of the tt'^, ry, and 
direct photon yields. 



simulated SVT, SSD, and IPC were underestimated by 
factors of K = 2, 2, and 1.2, respectively, compared to 
that in the real data [IS]. In simulations, the photon 
conversion probability in these detectors was in the range 
^conv = 0.3-3.3%. To account for the missing material 
in the GEANT model, the photon spectra were corrected 
by factors c/^^^ = (1 - E Pconv)/ (1 - E '«^'conv), with the 
values of 1.06 ± 0.02 and 1.03 ± 0.02 for the p + p and 
d + Au data, respectively. Using Eq. ([s]), this corresponds 
to correction factors of cj,gg = 1.12 ± 0.03 {p + p) and 
1.07 ± 0.03 {d + Au) for the 7r° spectra. Because the 
photon attenuation length in most absorbers rapidly ap- 
proaches a constant for energies larger than «100 MeV, 
the correction factors were assumed to be independent of 
the photon p^- 



IV. NEUTRAL PION AND ETA MESON 
ANALYSIS 

The 7r° and 77 were identified by their decays 

7]-'' — >. and 77 — > 77. 

These decay modes have branching ratios of 0.988 and 
0.392, respectively [50]. The 7r° lifetime is r = 8.4 x 
10"^'' s, corresponding to a decay length cr — 0.025 ^m. 
The 77 lifetime is even shorter (7 x 10~^^ s). Therefore, 
we assumed that the decay photons originated from the 
primary vertex. For each event, the invariant mass 



H. Photon conversions 

A separate study was done to determine the degree to 
which the GEANT geometry described the distribution of 
material in the real STAR detector, and the correspond- 
ing correction factors cioss were extracted to account for 
any differences. 

The photon conversion probability Pconv a-s a function 
of the depth d traversed in a material is given by 



-Pconv = 1 - cxp {-d/do) 



(4) 



where do is the mean free path of the photon in that 
material. The probablity that a 7r° was not detected 
because at least one of its decay photons has converted is 



-Ploss 2Pconv(l -Pconv) ~^ 



(5) 



The tt" losses due to conversions were in principle taken 
into account in the simulations mentioned in the sections 
below, because the material traversed by the photons was 
included in the GEANT model of the detector. However, 
it was observed that the simulation failed to reproduce 
the number of photon conversions in the inner tracking 
system (SVT, SSD), and in the TPC Inner Field Cage 
(IFC), all of which have a very complicated geometry of 
silicon sensors, readout electronics, and support struc- 
tures lini mi HE]- The number of conversions in the 



= ^2EiE2{l-cosi^) 



(6) 



was calculated for all pairs of photons detected in the 
BEMC. Here Ei and E2 are the energies of the decay pho- 
tons and ij} is the opening angle in the laboratory system. 
The reconstructed masses were accumulated in invariant 
mass spectra, where the tt*^ and the 77 showed up as peaks 
around their nominal masses (777,^0 = 0.135 GeV/c^ and 
772^ = 0.547 GeV/c^). These peaks were superimposed on 
a broad distribution of combinatorial background, which 
originated from photon pairs that were not produced by 
the decay of a single parent particle. 



A. Asymmetry of photon pairs 

The energy asymmetry of the two-body decay of neutral 
mesons is defined as 



Z = 



\E^-E-2\ 
El + E2 



(7) 



From the decay kinematics it follows that the probability 
for a given Z~^^ is independent on Z^-y. Figure [5] shows 
the distribution of the asymmetry of photon pairs recon- 
structed \n p -\- p data, including both tt^ and 77 signals 
and background. In the MinBias data the distribution is 
not flat because of the acceptance effects — photons from 
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FIG. 5: (color online) The energy asymmetry Z-,^ of photon 
pairs reconstructed in p + p data (symbols) and in Monte 
Carlo simulation (histograms) for various triggers, normalized 
to unity for each trigger. 
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FIG. 6: (color online) Distribution of vr" candidates obtained 
from the p + p data, as a function of (a) rj and (b) ip, and (c) 
Pt and (d) m^j. 



an asymmetric decay have a large opening angle and one 
of them is likely to escape the BEMC. It is also seen that 
the HighTower energy threshold biases the asymmetry to- 
wards larger values, because it is easier for an asymmetric 
decay to pass the trigger. The corresponding asymmetry 
distributions obtained from the Monte Carlo simulation, 
which represented the pure 7r° signal, are also shown. The 
details of the simulation are given in section |IV G| below. 
Asymmetries observed in the simulation are in general 
agreement with those in the real data, considering the 
presence of background in the data. 

In this analysis, the tt^ and rj candidates were only 
accepted if the asymmetry was less than 0.7. This cut 
rejected very asymmetric decays, where one of the BEMC 
points had low energy. It also rejected a significant part 
of the low-mass background (this background will be 
described in section [iV D ) . The asymmetry cut improved 
r° signal-to-background ratio by a factor of «1.5. 



the 



B. Kinematic distributions 

For each tt^ candidate, the pseudorapidity 77, the az- 
imuth (p, the transverse momentum px, and the invariant 
mass m^^ [Eq. ([61] were calculated. Figure [6] shows the ry, 
ip, Pt, and m-y.y aistributions of the tt^ candidates in the 
p+p data. For the d + Au data these distributions are 
similar to those shown for p + p. The corresponding 77, 
(fi, and Pt distributions of the tt'^'s reconstructed in the 
simulation are also shown. 

The r] distribution shows the decrease of the calorimeter 
acceptance at 77 = and rj = 1, because it is likely that 
one of the decay photons at the calorimeter edges escapes 
detection. The asymmetry of the 77 distribution is due 
to the fact that the calorimeter half-barrel is positioned 
asymmetrically with respect to the interaction region. 
The azimuthal dependence of the calorimeter acceptance 



was caused by failing SMD modules (the data used in 
this paper are from the early years of detector operation, 
in which such failures occured relatively frequently), as 
well as by dead and hot towers. The gross features of the 
data reflecting the calorimeter acceptance are reasonably 
well reproduced by the pure tt" Monte Carlo simulations. 

Figure ^c) shows the pt distribution of the pho- 
ton pairs, separately for the MinBias and for the two 
HighTower datasets. The HighTower trigger threshold 
effects are reasonably well reproduced in the simulation. 
It is seen that using the HighTower triggers significantly 
increased the rate of tt" candidates at high pT- 

The px-integrated invariant mass distribution in 
Fig. clearly shows the tt*^ and 77 peaks, superimposed 
on a broad background distribution. This background 
has a combinatorial and a low-mass component, discussed 
in detail in the two following sections. 



C. Combinatorial background 

The combinatorial background in the invariant mass 
distribution originated from pairs of photon clusters that 
were not produced by a decay of a single particle. To 
describe the shape of the combinatorial background, we 
used the event mixing technique, where photon candi- 
dates from two different events were combined. To avoid 
the mixing of different event topologies, the data were 
subdivided into mixing classes based on the vertex po- 
sition, BEMC multiplicity, and trigger type (MinBias, 
HighTower-1, and HighTower-2). 

Figure [7|[ a) shows, as an example, the invariant mass 
distribution in the 4 < pt < 5 GeV/c bin, obtained from 
the p + p HighTower-1 data, together with the combina- 
torial background obtained from the event mixing. The 
mixed-event background distribution was normalized to 
the same-event distribution in the invariant mass region 
0.9 < m^^ < 1.2 GeV/c2. 
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FIG. 7: The same-event invariant mass distribution (crosses) 
and combinatorial background (histogram) observed in one pr 
bin of the p + p data, (a) Background estimated from random 
event mixing and (b) from a linear combination of random and 
jet-aligned mixing, and (c) background-subtracted distribution. 
The shaded areas indicate the regions where the mixed-event 
background was normalized to the same-event distributions. 



There is still some residual background in the interval 
0.2 < m-y-y < 0.4GeV/c^. This background is due to 
correlation structures (jet structures) in the event, which 
are not present in the sample of mixed events. In order 
to preserve jet-induced correlations, the jet axes in both 
events were aligned before mixing [38j , as described below. 

To determine the rj-if position of the most energetic 
jet in every event, the cone algorithm was used [51j . The 
mixed-event tt'^ candidates were constructed by taking two 
photons from different events, where one of the events was 
displaced in i] and ip by A77 = ry2 — ??i and Aip = (p2 — ipi, 
respectively. Here 771 2 and (pi^2 are the jet orientations 
in the two events. 

Figure |8] shows a schematic view of two superimposed 
events where the jet axes are aligned. In order to minimize 
acceptance distortions, the events were divided into ten 
mixing classes in the jet 77 coordinate. By mixing only 
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FIG. 8: A schematic view of two superimposed events, where 
the jet axes are aligned to preserve the jet-induced correlations 
in the mixed event. 



events in the same class, the shift A77 was limited to 0.1. 
Because the calorimeter has a cylindrical shape, the shift 
in ip did not induce any significant acceptance distortion. 

A side effect of this procedure was that correlations 
were induced if there was no real jet structure, because 
the jet finding algorithm then simply picked the most 
energetic track in the event. To reduce possible bias intro- 
duced by such correlations, the combinatorial background 
was taken to be fully random for < 1.2 GeV/c and fully 
jet-aligned for px > 10 GeV/c. Between these values, the 
random component decreased linearly with increasing px- 
We assigned a systematic uncertainty of 10% to the ran- 
dom background fraction, which resulted in a systematic 
uncertainty of 5% of the 7r° and 3.5% of the t] yields. 

Figure [rj^b) shows the same invariant mass spectrum as 
that shown in panel (a), with the background estimated 
by the combined random and jet-aligned event mixing. 
The mixed-event background was normalized to the same- 
event distribution in the ranges 0.3 < m^-,, < 0.4 and 
0.9 < m^j < 1.2 GeV/c^. By changing the subtracted 
background within its normalization uncertainty, we ob- 
tained another component of a systematic error of the tt" 
and T] yields, which was found to increase with increasing 
Pt from 0.5 to 3% for the 7r° and from 10 to 50% for the rj. 

Figure [tJ^c) shows the background-subtracted distribu- 
tion. It is seen that there is still a residual background 
component at invariant mass m^^ < 0.1 GeV/c^. The 
origin of this background is described in the next section. 



D. Low-mass background 

Random fluctuations in the SMD signals occasionally 
generate a double-peaked hit structure, in which case the 
clustering algorithm incorrectly splits the cluster. These 
random fluctuations enhance the yield of pairs with mini- 
mal angular separation and thus contribute to the lowest 
di-photon invariant mass region, as can be seen in Fig.jTf^c). 
However, at a given small opening angle, the invariant 
mass increases with increasing energy of the parent parti- 
cle, so that the low-mass background distribution extends 
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FIG. 9: The simulated low-mass background distributions 
from erroneous splitting of single photons, in three bins of 
the reconstructed pair px- The distributions extend to larger 
invariant masses with increasing pr and move into the w'^ 
region (shown for pr ~ 15 GeV/c). 
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to larger values of with increasing photon pT- 

The shape of the low-mass background was obtained 
from a simulation as follows. Single photons were gen- 
erated with flat distributions in — tt < Lp < -I- tt, 
-0.2 < ri < 1.2, and < pr < 25 GeV/c. These pho- 
tons were tracked through a detailed description of the 
STAR geometry with the GEANT program. A detailed 
simulation of the electromagnetic shower development in 
the calorimeter was used to generate realistic signals in 
the towers and in the SMD. The simulated signals were 
processed by the same reconstruction chain as the real 
data. Photons with more than one reconstructed cluster 
were observed, and m^^ and pt of such cluster pairs were 
calculated. The m^^ histograms were accumulated, with 
each entry weighted by the pr spectrum of photons in the 
real data, corrected for the photon detection efficiency. 

Figure |9] shows the low-mass background distributions 
in three bins of the reconstructed pair px- It is seen that 
the distributions indeed move to larger invariant masses 
with increasing pT and extend far into the tt*^ window at 
high Pt- For this reason, it was not possible to estimate 
the amount of this background from a phenomenological 
fit to the data, and we had to rely on the Monte Carlo 
simulation to subtract this background component. 

The second significant source of BEMC clusters that 
passed the CPV cut was the neutral hadrons produced 
in the collisions, mostly antineutrons above 2 GeV/c. To 
account for the additional low-mass background from 
these hadrons, simulations of antineutrons were performed 
in the same way as of photons, and the reconstructed 
invariant mass distribution was added according to the 
realistic proportion n/7inch where the antineutron yield 
was estimated as described in section IVB] below. 

The low-mass background was normalized by matching 
the observed pt spectrum of the clusters between sim- 
ulation and data. This removal procedure worked well, 



FIG. 10: The invariant mass distribution observed in one pr 
bin of the p + p data (a) before and (b) after the low-mass 
background subtraction. 



and Fig. |10| shows the invariant mass spectra and the low- 
mass background component, and the final background- 
subtracted spectrum for the p-\-p HighTower-I data. The 
normalization uncertainty of the low-mass background 
contributes to the systematic uncertainty of the tt*^ cross 
section, and reaches 15% at the high-px end of the spec- 
trum. 



E. Peak position and width 



Figure 11 'a) shows the background-subtracted 



distribution in the region 4 < pt^ < 5 GeV/c obtained 
from the p + p HighTower-1 data (symbols) , together 
with the corresponding distribution from the detector 
simulation (histogram). In order to compare the real and 
simulated m^^ distributions for all bins in p^ and for all 
datasets, we estimated the position and width of the peaks 
using Gaussian fits in the peak region. Figure [Tl|b) shows 
the peak positions obtained from the fit to the p + p data. 
It is seen that the peak position shifts towards higher 
masses with increasing pT- This shift is a manifestation 
of the bin migration effect that originates from statistical 
fluctuations in the calorimeter response. Due to the 
steeply falling pT spectrum, the energy resolution causes 
a net migration towards larger px- Since larger values of 
Pt imply larger values of m^^, the migration effect biases 
the invariant mass peak towards larger values. 

An additional peak shift at the largest values of p-r 
is caused by the SMD strip granularity, which imposes 
a lower limit on the opening angle of the reconstructed 
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FIG. 11: (color online) (a) Invariant mass spectrum of 7r'''s 
reconstructed in the simulation (histogram) in comparison to 
the p + p HighTower-1 data (symbols) in 4 < pT < 5 GeV/c 
bin. (b) Peak position and (c) width in the real data (solid 
symbols) and in the simulation (open symbols). Horizontal 
line in panel (b) shows the true tt'' mass. 



FIG. 12: (color online) (a) Invariant mass spectrum of 77 
mesons reconstructed in the simulation (histogram) in com- 
parison to the p + p HighTower-1 data (symbols) in 4 < pr < 
5 GeV/c bin. (b) Peak position and (c) width in the real 
data (solid symbols) and in the simulation (open symbols). 
Horizontal line in panel (b) shows the true rj mass. 



photon pairs. The minimal SMD cluster separation in 
each dimension that can be resolved by the cluster finder 
is 1.5 strips, and most clusters contain at least two or three 
strips. Therefore, the pair reconstruction is less efficient 
for the symmetric decays with the smallest opening angles 
at pt ^ 10 GeV/c. This leads to an increased average 
opening angle and m^^ of the reconstructed photon pairs 
from tt" decays. 

The peak position observed in the data is larger than 
that found in the simulations by (3.5 ± 0.6)%, on average. 
This difference could be caused by the global energy 
scale of the BEMC towers being off by a similar amount. 
We already accounted for this possibility by assigning a 
systematic uncertainty of 5% to the BEMC calibration 



constants (see section HIE). 

c) compares the n'^ peak width in the data 



Figure 



II 



and in the simulation, and it is seen that the peak width in 
the data is larger than that in the simulation by (25±2)%, 
on average. This is a sufficiently good agreement for this 
analysis, because the 7r° and r] yields were counted in the 
mass windows that were adjusted in each px bin to cover 
the entire signal peak. 

The peak shape of the rj meson, as well as its position 
and width, are shown in Fig. [T2| as a function of the 
reconstructed pt- The peak position in the data is larger 
than in the simulations by (5.1 ± 1.2)%, and the width 
by (34 ± 11)%, on average. The difference in the peak 
position is similar to the ir'^ case above, which supports 
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the possibility of both being caused by a small systematic 
offset in the BEMC calibration. The observed level of 
agreement between data and simulations is considered to 
be sufficient for this analysis. 



F. Invariant yield extraction 

The invariant yield of the tt" and 77 mesons per MinBias 
collision, as a function of px, is given by 



E 



dp-^ Pt dpT dy dtp 2ttpt dpT dy 



(8) 



where in the last equality integration over the full 27r 
azimuthal coverage of the STAR detector is performed. 
Using the experimentally measured quantities, the invari- 
ant yield was calculated as 



E 



(9) 



where: 



• y is the raw yield measured in the bin centered at 
Pt and y; 

• ApT is the width of the px bin for which the yield 
was calculated; 

• Ay is the rapidity range of the measurement; in this 
analysis Ay = 1 for all data points, except for the 
T] yields at < 3 GeV/c, where the correction for 
the difference between rapidity and pseudorapidity 
reached 7%; 

• A'tiig is the number of triggers recorded; 

• -ftTtrig is the trigger scale factor; -fCtrig = 1 for the 
MinBias events and >1 for the HighTower events; 
the product -/Vtrig-ftTtj-ig is the equivalent number of 
MinBias events that produced the yield Y; 

• Evert is the vertex finding efficiency in MinBias 
events; 

• cJJ^gg is the correction for the missing material in the 
simulation; 

• e^cc is the BEMC acceptance and efficiency correc- 
tion factor; 

• ejpv is a correction for random TPC vetoes; 

• — T-yy/r is the branching ratio of the di-photon 
decay channel (0.988 for tt" and 0.392 for rj [50 ). 

The raw n'^ and 77 yields were counted in the px- 
dependent m-^-y windows that contained the peaks. The 
low-mass border of the tt" peak region was taken to be a 
linear function ofpx, common for all datasets and triggers. 
This cut was optimized to capture most of the yield and as 



little of low-mass background as possible. The high-mass 
border also linearly increased with pT in order to cover 
the asymmetric right tail of the peak. Similarly, the rj 
peak region was a p^-dependent window that captured 
most of the signal. For completeness, we give below the 
parametrization of the tt^ and 77 windows: 



75 + l.lpT < m^^(7rO) < 250 + S.Spr MeV/c^ 



350 + 3.3pT < niy^iTj) < 750 



MeV/c^ 



(10) 



where px is measured in GeV/c. The stability of the yields 
was determined by varying the vertex position cut, the 
energy asymmetry cut, and the yield integration windows. 
From the observed variations, a point-to-point systematic 
error of 5% was assigned to the 7r° and 77 yields. 

Within each trigger in the p + p data, the tt" signal 
significance decreased from sa34 to w6 standard deviations 
with increasing pT, because of the corresponding reduction 
in statistics. In the d + Au data, the same trends were 
observed, but the significance was lower than in the p + p 
data by a factor of 1.9, on average, which is mainly caused 
by the lower integrated nucleon-nucleon luminosity in 
these data. The significance of the f] signal was between 
18 and 2.5 standard deviations in the p + p data, and 
between 5.5 and 1.0 standard deviations in the d + Au 
data. 



G. Acceptance and efficiency correction 

To calculate the detector acceptance and reconstruction 
efficiency correction factor e^^.^., a Monte Carlo simula- 
tion of the detector was used. The 7r° decay photons 
were tracked through the STAR detector geometry using 
GEANT. The simulated signals were passed through the 
same analysis chain as the real data. 

The tt'^'s were generated in the pseudorapidity region 
— 0.3<r7< +1.3, which is sufficiently large to account 
for edge effects caused by the calorimeter acceptance 
limits of < 77 < 1. The azimuth was generated flat in 
— TT < (f < +77. The Pt distribution was taken to be 
uniform up to 25 GeV/c, which amply covers the measured 
pion Pt range of up to 17 GeV/c. The vertex distribution 
of the generated pious was taken to be Gaussian in z with 
a spread of cr = 60 cm and centered at z = 0. 

The generated 7r°'s were allowed to decay into two pho- 
tons, TT*' — > 77. The GEANT simulation accounted for all 
interactions of the decay photons with the detector, such 
as conversion into e'^e~ and showering in the calorimeter 
or in the material in front of it. 

To reproduce a realistic energy resolution of the 
calorimeter, an additional smearing had to be applied 
to the energy deposits calculated by GEANT. In all simu- 
lations, a spread of 5% was used to reproduce the p + p 
data and 10% for the d + Au data. 

To reproduce the pt spectrum of pions in the data, 
each Monte Carlo event was weighted by a pT-dependent 
function. This weighting technique allowed us to sample 
the entire pt range with good statistical power, while 
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FIG. 13: (color online) Acceptance and efficiency correction 
factor Eacc for (a) tt" and (b) i], calculated from the Monte 
Carlo simulation for the p + p data. 



reproducing the bin migration effect caused by the finite 
detector energy resolution. An NLO pQCD calcula- 
tion [32] provided the initial weight function, which was 
subsequently adjusted in an iterative procedure. 

The time dependence of the calorimeter acceptance 
during data taking was recorded in database tables that 
were used in the analysis. In order to reproduce this time 
dependence in the simulation, the generated events were 
assigned timestamps that followed the timeline of the 
real data taking. In this way, the geometrical acceptance 
of the calorimeter (mean fraction of good towers) was 
reproduced in the Monte Carlo with a precision of «0.5%. 

In the analysis of real data, we used vertices re- 
constructed from TPC tracks and those derived from 
BBC time-of-flight measurements. The former have sub- 
millimeter resolution, whereas the latter have a precision 
of only «40 cm. To account for the BBC vertex resolu- 
tion, 35% of the generated pions in the p + p MinBias 
data had their point of origin artificially smeared in the 
z direction. No such smearing was applied to the other 
simulated data, where BBC vertex was not used. 

The acceptance and efficiency correction factor was 
calculated from the simulation as the ratio of the raw tt" 
yield reconstructed in a px bin to the number of simu- 
lated pions with the true pr generated in that bin. This 
was done separately for each trigger, using the same tt" 
reconstruction cuts as was done in the real data analysis. 
In particular, the reconstructed value of pseudorapidity 
was required to fall in the range < 77 < 1 in both the 
data and the simulation, while in the latter the generated 
value of 77 was also required to fall in this range. As an 
example. Fig. 13 shows the 7r° and 77 correction factors 
for the three triggers in p + p data. 

The difference between the MinBias and HighTower 
correction factors was caused by the SMD requirement 
in the HighTower data, which was absent in the MinBias 
data. The absence of the SMD information reduced 
the 7r° reconstruction efficiency in the MinBias data at 
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FIG. 14: Acceptance and efficiency correction factor Eacc for 
the p+p HighTower-1 data, with standard set of cuts (solid 
symbols) and with SMD quality cut removed (open symbols). 



> 3 GeV/c, where the decay photons were separated 
by less than two towers. The 77 reconstruction is only 
affected by this at larger values of pT ■ 

The effect of the SMD quality requirement of having 
at least two adjacent strips in a cluster is illustrated in 
Fig. [T4j which shows the correction factor calculated for 
the p + p HighTower-1 data with (solid symbols) and 
without (open symbols) the SMD quality requirement. 
This requirement reduced the number of accepted tt" 
candidates by «45%. This explains the difference between 
the HighTower-1 and HighTower-2 (no SMD quality cut) 
correction factors at high pT seen in Fig. 
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The current simulation framework poorly reproduces 
the shower shapes in the SMD at the low incident photon 
energies. To account for residual bias after applying the 
SMD quality cut, we assigned a systematic uncertainty 
to the HighTower-1 cross section, which decreases from 
15% at PT = 4 GeV/c to zero at pt = 7 GeV/c. 

To determine a dependence of the acceptance correction 
on the track multiplicity M , and thus on the centrality, 
we analyzed a sample of generated tt^'s embedded into 
real d + Au events. No significant centrality dependence 
was found. Therefore, the same correction factors were 
applied to the different centrality classes in the d-|-Au data. 
The dependence of the efficiency on the locally higher 
multiplicity in jets was investigated in a pythia [52j 
simulation, and no significant difference in the efficiency 
was observed relative to a single particle simulation. 



H. HighTower trigger normalization 

We have shown in Fig. |6] the pt distribution of tt^ candi- 
dates for the p+p MinBias, HighTower-1, and HighTower-2 
data. To normalize the HighTower spectra to those of 
the MinBias, pT-independent scale factors were applied. 
These scale factors were estimated as the ratio of observed 
MinBias to HighTower event rates. 



(11) 



Here A^mb and A'ht are the numbers of MinBias and 
HighTower triggers that passed the selection cuts, Smb 
and are the online prescale factors adjusted on a run- 
by-run basis to accomodate the DAQ bandwidth, and 
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the sums are taken over all runs where both the MinBias 
and HighTower triggers were active. We obtained the 
values i^tiig = 4.67 x 10^ and 1.96 x 10** for the p + p 
HighTower-1 and HighTower-2 triggers, respectively, and 
Ktrig = 2.87 X 10^ and 2.86 x 10^ for the d + An triggers. 

To check the scale factors, the HighTower software 
filter, which simulated the hardware trigger, was applied 
to the MinBias data. The scale factors were obtained as 
the ratio of the total number of MinBias events to the 
number that passed the filter. To obtain a more precise 
HighTower-1/ HighTower-2 relative normalization factor, 
the software filter was applied to the HighTower-1 data. 
The results from the two methods agreed within 3% for 
HighTower-1 data and within 5% for HighTower-2 data. 
These numbers were taken as the systematic uncertainties 
of the HighTower normalization factors. 

The difference between vertex finding efficiencies in 
MinBias and HighTower data was effectively absorbed 
in the scale factor Knig. Therefore, the vertex finding 
efficiency correction was applied to the scaled HighTower 
data, as well as to the MinBias data. 



I. Bin centering scale factors 

To assign the yield measured in a pr bin to a single 
Pt value, the procedure from Ref. [S3] was applied. The 
variation of the yield within a bin was approximated by 
the function /(pr) = Aexp{—BpT)- The measured yield 
in the bin was assigned to the momentum p^ calculated 
from the equation 
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The procedure was repeated, taking p^ as the abscissa, un- 
til the p^ values were stable (typically after 3 iterations). 

To facilitate the comparison of results from the various 
datasets, the yields were scaled to the bin centers by the 
ratio iiTbin = /(Pt)//(pt), where pt is the center of the 
bin. The statistical and systematic errors were also scaled 
by the same factor. 

To estimate a systematic uncertainty introduced by this 
procedure, we changed the functional form of /(pt) either 
to a local power law in each bin or to a global power law 
in the full pt range. The observed variation in i^bin was 
below 1.5% in most pr bins. 



J. Fully corrected yields 

The fully corrected 7r° and f] invariant yields per 
MinBias event [Eq. ^] are shown in the top panels of 
Figs. [T5|and |16[ respectively. The curves represent power 
law fits to the data, of the form 
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FIG. 15: (color online) (a) Invariant yield of tt" per MinBias 
p + p and d + Au collision. Curves are the power law fits 
given in the text. Invariant yield divided by the fit to the (b) 
p + p and (c) d + Au data. The error bars are statistical and 
brackets in the lower panels are the systematic uncertainties. 



The ratios between the data and the fits are shown in 
the lower panels. The brackets in the lower panels are 
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the systematic uncertainties (see section IV K below), 
which are partially correlated between different trigger 
datasets. The agreement between the data taken with 
the different triggers is satisfactory, although a small 
systematic difference between the d + Au HighTower- 
1 and HighTower-2 yields (1.9 standard deviations, on 
average) was observed. 

For the calculation of the final cross sections and their 
ratios, the data from three triggers were merged. The 
HighTower-1 data were used in the MinBias-HighTower-1 
overlap bins, because the MinBias data in those bins al- 
most entirely represent a subset of HighTower-1 events, 
selected by an online prescale factor S'mb- Similarly, 
the HighTower-2 data were used in the HighTower-1- 
HighTower-2 overlap bins. 
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TABLE III: Systematic error contributions. The classifications A, B, C, and N are defined in the text. The error contributions 
to the 7r° cross section, the ?;/7r° ratio, Rcp, RdA, and the inclusive photon yield are indicated in the respective columns. 



Source 


Type 


Value at low (high) pr (%) 


E(f(j/dp^ 


r?/7r° Rcp 


RdA 


7inc 


Combinatorial background 


A 

A 


U.O (o) 


+ 


+ + 


+ 




Mixed- event background 




p: 



1 

+ 


1 

+ 






Low- mass background 


\^ 


1 f'\ K\ 

1 (15) 


1 

+ 


+ 






Random vetoes 


1S.T 

IN 




+ 


+ 


+ 


+ 


HighTower normalization 


a 


3 (5) 


+ 


+ 


+ 


+ 


/inaiybio culo 


A 




1 




1 




Conversion correction 


B 


3 for tt'' and 77, 2 for 7inci 


+ 




+ 


+ 


Tower energy scale 


B 


15 (35) in p + p, 10 (40) in d + Au 


+ 




+ 


+ 


SMD simulation 


C 


15 (0) at pT = 4 (7) GeV/c 


+ 








SMD energy scale 


B 


3.5 (1) at pT = 4 (11) GeV/c 


+ 




+ 




Bin centering 


C 


1.5 


+ 








Vertex finding efficiency 


N 


1 in d + Au MinBias 


+ 


+ 


+ 


+ 


MinBias cross section 


N 


11.5 in p + p, 5.2 in d + Au 


+ 




+ 


+ 


Glauber model (A'^coii) 


N 


5.3 in d + Au MinBias, 10.5 in Rcp 




+ 


+ 





K. Summary of systematic uncertainties 

The uncertainty of the calorimeter tower calibration 
was the dominant source of systematic uncertainty in this 
analysis. The uncertainty of the uncorrected yield Y{pt) 
was estimated from 



5Y{pt) 



dY{pT) 



dpi 



Spj 



(14) 



where Spr was taken to be 5% in the d + Au and p + p 
data (as derived from the electron calibration, see sec- 
tion HIE). This pT-dependent systematic uncertainty 
was, on average, 25% in both p + p and d + Au data. 

Another strongly py-dependent term in Eq. ^ is the 
acceptance and efficiency correction factor £acc, obtained 
from the Monte Carlo simulation. However, this term 
contributed much less to the cross section uncertainty, 
because the simulations used the inverse of real calibra- 
tion constants to convert the GEANT energy deposit in 
each tower to the ADC value. Therefore, the energies 
later reconstructed from the ADC values with the same 
constants were not sensitive to their fluctuations. This 
factor was sensitive to the tower calibration only through 
the HighTower threshold values, and contributed a sys- 
tematic uncertainty of 5% in the p + p data and 8% in the 
d + Au data, on average, correlated with the uncorrected 
yield uncertainty estimated above. 

The uncertainty of the SMD energy scale entered the 
analysis mainly due to threshold effects in the clustering, 
where the loss of a soft photon from the asymmetrically 
decayed ir'^ may change the reconstruction efficiency. This 
analysis does not depend on the absolute energy calibra- 
tion of the SMD, because the main energy measurement 
was obtained from the towers. Instead, we estimated the 
possible disagreement between the SMD scale in the data 
and in Monte Carlo simulation to be below 20%. The 
resulting variation of £acc from this source was 3.5% in 



the HighTower data at the pion px — 4: GeV/c and less 
than 1% above 11 GeV/c. 

All systematic error contributions are summarized in 
Table [Till classified into the following categories: 



A point-by-point systematic uncertainty; 

B pT-correlated systematic uncertainty, but uncorrelated 
between p + p and d + Au datasets; 

C pT-correlated systematic uncertainty, also correlated 
between p + p and d + Au datasets; 

N normalization uncertainty, uncorrelated between p + p 
and d + Au datasets. 

Table |ni| also lists which measurement is affected by a 
given source of systematic error. 



V. DIRECT PHOTON ANALYSIS 

The traditional approach to measuring direct photon 
production in hadronic collisions uses isolation criteria. 
Photons from decays of highly energetic hadrons should 
be accompanied by other jet fragments. Therefore, one 
can reject those by requiring less than a certain amount 
of background energy in a cone around a photon can- 
didate [55]. However, prompt photon production be- 
yond leading order in pQCD cannot be separated un- 
ambiguously from photons from fragmentation processes, 
although the framework for applying isolation cuts in 
pQCD calculations is established and theoretical inter- 
pretation of experimental results is possible. In addition, 
the use of isolation cuts in the high-multiplicity environ- 
ment of heavy-ion collisions is not straightforward. 

As the present analysis is intended to provide a baseline 
measurement for heavy-ion collisions, we have chosen to 
use the method of statistical subtraction to obtain direct 
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FIG. 16: (color online) (a) Invariant yield of per MinBias 
p + p and d + Au collision. Curves are the power law fits 
given in the text. Invariant yield divided by the fit to the (b) 
p + p and (c) d + Au data. The error bars are statistical and 
brackets in the lower panels are the systematic uncertainties. 



photon yields. For this method, one obtains inclusive pho- 
ton spectra, which, in addition to the direct contribution, 
contain a large background of decay photons, doniinantly 
from tt" decays. An accurate measurement of tt'^ and 
heavier hadrons provides the necessary input to subtract 
the decay background. This method has been successfully 
used in heavy- ion reactions [561 157] . however, it does not 
provide event-by-event direct photon identification. 

The sample of photon candidates served as the main 
input to the direct photon analysis, as in the case of 
the reconstruction of the tt'^ spectrum, described in the 
previous sections. After subtracting the contamination by 
charged particles and neutral hadrons, the raw inclusive 
photon sample was corrected to account for the limited 
acceptance and the finite detector resolution. In parallel, 
the total yield of photons from 7r°, rj, and w(782) decays 
was simulated, assuming a phenomenological scaling law 
{rriT scaling) for the and a;(782) spectra. 

To exploit the fact that the inclusive photon and decay 
photon yields have many correlated uncertainties, we 



7incl 



Tdccay /'^^ 



(15) 



where the numerator equals the point-to-point ratio of 
the measured spectra of inclusive photons and 7r°'s, as 
a function oi pr, and the denominator is the simulated 
background contribution from decay processes divided by 
the parametrized 7r° yield. It follows that 



1 



7dii 



Tdccay 



(16) 



which serves as an indicator of a direct photon signal 
7dir {Rj > !)■ The absolute direct photon yields can 
subsequently be determined as 



7dir 



1-r: 



7incl, 



(17) 



where the systematic uncertainties, which canceled in the 
double ratio, have to be included again. 



A. Inclusive photons 

The reconstruction of the inclusive photon spectrum 
was in many ways similar to that of the tt". The un- 
corrected photon spectrum was extracted from the same 
data sample using identical event and photon candidate 
cuts. However, in the tt*^ analysis there were no rigorous 
constraints on the purity of the photon candidates, be- 
cause remaining contributions from charged particles and 
neutral hadrons were identified afterwards as the combina- 
torial background in the mass distributions. In contrast, 
the uncorrected inclusive photon yield linci was obtained 
from an explicit subtraction of such backgrounds, 



cand 1 



(18) 



where the correction terms Cq and C'± represent the frac- 
tional contamination by neutral hadrons and charged 
particles, respectively. The charged particle contamina- 
tion C'± was estimated in section III G and found to be 
smaller than 5%. 

The invariant yield of inclusivephotons was calculated, 
similarly to that of 7r°'s in Eq. (|9|), as 



E- 



incl 



dp3 2TrpTNtTigKtrig AprAy elccEi 



(19) 



Here e^cc is the single photon acceptance and efficiency 



correction factor, discussed in section V C below 



B. Neutral hadron background 



The term Co in Eq. ( 18 ) was defined as the number 



of reconstructed showers generated by neutral hadrons 
relative to the total number of showers in the photon 
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FIG. 17: (color online) Relative neutral hadron contribution 
to the photon candidate yields, Co- Filled area is the upper 
limit for the p + p HighTower-2 data sample, dashed curve 
is the upper limit in the extreme scenario that all photon 
candidates were a result of showering n's. 



Two natural limits to the contamination were consid- 
ered: (i) the hard upper limit Co < 1, which is not 
realistic, as it implies that the inclusive photon yield is 
zero, and (ii) the limit implied by the assumption that 
the direct photon signal is zero around the annihilation 
peak («2 GeV/c). In both cases, a scaling factor for Co 
was derived and subsequently applied to the HighTower-2 
data, as shown in Fig. |17[ For the further analysis, we 
have chosen the second estimate, calculated assuming 
that only background photons were detected in the range 
1 < Pt < i GeV/c, as the upper limit of the neutral 
hadron contamination Cq. This upper limit was found 
to be negligible in the range of the present direct photon 
measurement, px > 6 GeV/c. 

The n and interactions with the BEMC resulted 
in the smaller contamination than that of the n's, at all 
values of px ■ 



candidate sample. The STAR detector has no means of 
directly identifying neutrons and antineutrons. Therefore, 
this contamination was simulated using the measured 
(anti)proton spectra as input. 

The largest source of neutral contamination was the 
n annihilation in the calorimeter, for example n + p — t- 
27r"'"7r~7r'^. This initiates a shower that does not necessar- 
ily develop in the incident direction of n. Moreover, the 
available energy for the reaction products includes twice 
the rest mass of a nucleon (ss2 GeV). 

STAR has measured the p and p production in p + p 
and d + Au collisions [21] ■ The reported yields, however, 
were not corrected for the A and A feed-down, which 
is expected to have a contribution of Sa « 20% [23] , 
Therefore, the n yield was estimated as 

Yin) = (1 - Sa) Y{p) + ^A ^|^^^J°j y(p), (20) 

where the branching ratios are B{K — > titt") = 0.358 and 
B{K PTT-) = 0.639 HO]. 

To study the contamination of the photon candidate 
spectrum, approximately 3 x 10^ n's were generated with 
an exponentially falling px spectrum. This provided suf- 
ficient statistics at low px, where the n's constituted 
a significant source of contamination. The fluka pro- 
gram j58| was used to describe the particle transport and 
the interactions in the detector material. The parametriza- 
tions of the p and p yields were not only used to assign a 
weight to the Monte Carlo events, but also to determine 
the absolute contamination of the photon sample. The 
latter was divided by the number of reconstructed photon 
candidates to calculate the term Cq. 

The final contamination factor Co in the p + p data is 
shown in Fig. [17] In a limited px range, Cq appears to be 
larger than unity, which is not possible unless the asso- 
ciated systematic uncertainties are extremely large. The 
violation of this bound, as well as observed discrepancies 
between the three triggers, indicate the large uncertainty 
of the n simulation. 



C. Photon reconstruction efflciency 

We have calculated the photon acceptance and efficiency 
correction factor e^cc separately for events containing 7r° 
decay photons and for events containing only a single 
photon. The latter factor was applied to the fraction of 
the photon yield from all sources other than the tt" decay. 

To determine the acceptance correction eJcc ^ot the 
tt'^ decay photons, we used a GEANT-based Monte Carlo 
simulation of the STAR detector. The Monte Carlo 
events were weighted in such a way that the measured 
TT^ yield was reproduced. This is important because the 
photon acceptance depends on the degree of cluster merg- 
ing. This, in turn, depends on the opening angle of the 
decay photons, and thus on the momentum of the parent 
tt". Furthermore, the simulation included all the possible 
losses of photon candidates listed in section |IV[ except 
those associated with the invariant mass window and with 
the cut on the energy asymmetry Z-^^. One important 
effect affects the showers initiated by daughters of a high- 
px 7r°. Because there was no requirement on the relation 
between the reconstructed px and the Monte Carlo input 
Px, the correction implicitly accounted for events in which 
one of the two decay photons remained unidentified and 
the total energy was assigned to a single cluster. Such 
merging of photon showers constituted the main differ- 
ence between the reconstruction efficiencies of tt^ decay 
products and single photons. 

Similarly to the above, the e^cc factor for single photons 
was determined using a Monte Carlo sample of wl x 10^ 
events. Each event contained a single photon, uniformly 
distributed in azimuthal angle —TT<(p< + ir, pseudo- 
rapidity —0.3 < 77 < -|- 1.2, and transverse momentum 
< < 20 GeV/c. Events were weighted with a function 
determined from the spectrum of photons from decaying 
hadrons other than the 7r°, as well as from that of the 
direct photons. However, we will demonstrate in sec- 
tion jV E| that the shape of the decay photon spectrum 
and the measured tt" spectrum were very similar, at least 
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TABLE IV: Dominant hadronic decay contributions to the 
inclusive photon yield [S^ . 



Decay 


Branching ratio (%) 


tt" — > 77 


98.80 


tt" — > e^e~7 


1.20 


— )■ 77 


39.23 


77 — > 7r^7r^7 


4.78 


rj — > e^e~7 


0.49 


a;(782) ^ TV°-y 


8.69 



ues considered in this analysis, the decay photon efficiency 
even exceeded unity. 



FIG. 18: (color online) (a) Acceptance and efficiency factor 
£acc for photons exclusively from the decay tt" — >■ 77. The 
observed rise of the efficiency at high pr is caused by the 
merging of the tt" decay daughters, (b) The eZcc factor for 
single photons, used to correct the fraction of the photon 
sample that exceeded the simulated tt" decay contribution. 



for the Pt range of this analysis. Although the direct pho- 
ton spectrum was expected to exhibit a slightly different 
Pt dependence, varying the input spectrum correspond- 
ingly did not yield quantitatively different results. 

Finally, we have implemented a correction to the mea- 
sured photon yields as follows. The yield of photons that 
originated from the decay tt'^ — > 77 was determined from 
the measured tt" spectrum. This part of l^nci was cor- 
rected with the e2cc factor calculated for the 7r° decay 
photons. The remaining part of i^nci was assumed to con- 
sist of single photons that were not correlated with the 
other photon candidates in the event, and was corrected 
with the single-photon eJcc factor. This assumption was 
based on the observation that the reconstruction efficiency 
for photons from the decay — >■ 77, which is the second 
largest source of decay photons («15%), was equal to that 
of single photons, because the opening angle between the 
two daughter photons is large enough that both are never 
incident on the same calorimeter tower. 

Figure [18] shows the acceptance and reconstruction effi- 
ciency ej^^ for the 7r° decay photons and for single photons 
for the p + p MinBias, HighTower-1, and HighTower-2 
data. The two HighTower results were found to be very 
similar in the low-pT region, where the angular separa- 
tion of the decay photons was still large, compared to 
the size of a BEMC tower. However, at higher pT the 
two photons are difficult to separate, particularly in case 
of the most symmetric decays. When two photons were 
merged, the remaining photon candidate was erroneously 
assigned the energy sum of both showers. This led to 
significantly larger reconstruction efficiency, compared to 
that for single photons. Eventually, at the largest pT val- 



D. Fully corrected inclusive yields 

The inclusive photon yield llnci was derived by sub- 
tracting the charged and neutral backgrounds from the 
yield of raw photon candidates [Eq. (18)]. The contami- 
nation by charged particles was subtracted according to 
the procedure explained in section [ill G[ but the neutral 
hadron correction proved to be difficult. Although an 
upper limit for the contamination fraction Cq was derived 
in section [V B| we did not find any means to reduce the 
associated systematic uncertainty on Cq to a level where 
a meaningful subtraction could be performed for MinBias 
data. In the HighTower data, the upper limit on the 
contamination fraction vanishes at higher values of pt- 
Therefore, our final results were obtained in the range 
6 < Pt < 15 GeV/c, and the photon candidates obtained 
from the MinBias data were discarded. 

Figure [T9| shows the corrected inclusive photon spectra 
in p + p and d + Au collisions without the subtraction 
of the neutral hadron contribution. The lower panels 
show the data divided by the corresponding power law 
fits. A small systematic difference between the spectra 
from d-l- Au HighTower-1 and HighTower-2 collisions (1.2 
standard deviations, on average) was observed. However, 
the measured HighTower-2 yields is statistically more 
significant, because the HighTower-1 events that contained 
photons in the range 6 < pt < 10 GeV/c were a subset 
of the HighTower-2 data. We reconstructed 17684 (3738) 
photon candidates from the HighTower-2 (HighTower-1) 
d-|-Au data in that pt range. The final direct photon cross 



sections presented in section VI were obtained exclusively 
from the HighTower-2 data. 



E. Background from hadronic decays 

The photon yield from hadronic decays was determined 
with a simulation of the decay processes listed in Table |IV] 
The other possible contributions, from processes such as 
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FIG. 19: (color online) (a) Inclusive photon invariant yield 
per MinBias p + p and d + Au collision. Curves are the power 
law hts given in the text. Invariant yield divided by the ht 
to the (b) P + P and (c) d + Au data. The neutral hadron 
contamination was not subtracted (see text) and is expected 
to be significant at 2 < pt ^ 4 GeV/c. The error bars are 
statistical and brackets in the lower panels are the systematic 
uncertainties. Data points for the overlapping pr bins in the 
lower panels are horizontally displaced for clarity. 



r]' -+ /9°7, were found to be negligible (<1%). A fit of the 
measured 7r° yield in the range 4 < pT < 15 GeV/c to the 
form ~(1 +pt)~" served as an input to the simulation. 
The fit yielded a = 9.1 ± 0.1 and 9.0 ± 0.1 for d + Au 
and p + p collisions, respectively. The normalization is 
irrelevant because it cancels in the ratio 7dccay/7'''^. 

To estimate the yields of the r] and w(782), we used the 
fact that these scale with the 7r° yields when expr essed in 
terms of the transverse mass ttit = y/rn? + p"^ instead 



of pt [13 . For the 77 spectra, we used the scaling 
ratios R^/^ = 0.46 ± 0.05 for p + p and 0.44 ± 0.08 for 
d + Au data, as followed from our measurement of i] 
production. In case of tj(782), we used R^/tt = 1-0 ± 0.2, 
in agreement with recent measurements at RHIC |62| . 
The estimated 77 and a;(782) yields relative to the 7r° yield 
are shown in Fig. [20] 

shows the ratio 7dccay for d-l-Au collisions. 



Figure 
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FIG. 20: The estimated yield of rj and a;(782) mesons in d-f Au 
collisions, relative to the measured tt" yield, determined from 
the rriT scaling as described in the text. 
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The curves represent the contributions of the tt , 77, and 



FIG. 21: The simulated number of photons per input pion 
7dGcay/7i"'' from hadronic decays in d -I- Au collisions, as a 
function of pr- The included decay processes are listed in 
Table IIVI The vertical line indicates the lower limit of the 
Pt range that was used to fit the tt" spectrum. The value 
of 7dccay/7r° below Pt = 4 GeV/c is, therefore, less accurate, 
however, it was not used for further analysis (see text). 



w(782), and tlie total decay photon yield, each divided by 
the parametrization of the measured tt*^ spectrum. The 
normalization uncertainty cancels upon taking this ratio. 
The uncertainty due to the shape of the tt'^ spectrum and 
to the TTiT scaling factors was estimated by varying the 
fitted exponents and the scale factors by their errors. 



F. Summary of systematic uncertainties 



All systematic error contributions to the double ratio 
R^ [Eq ([T5|)] are summarized in Table 



V[ Expressing 

the direct photon yield in terms of a double ratio gives a 
large reduction in the systematic error, since the contri- 
bution from the BEMC energy scale uncertainty cancels. 
Consequently, the largest sources of uncertainty are those 
associated with the tt" yield extraction and with the SMD 
energy scale. The latter leads to the i?^ variation of 12%, 
independent of pt for px > 6 GeV/c. 

The beam background observed in the d + Au data 
has a larger effect on single-photon analysis than on the 
tt'^ reconstruction, since the background-induced showers 
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TABLE V: Systematic error contributions for the double ratio 
Rj. The classifications A and B are defined in section [V] 



Source 



Type Value at low (high) pr (%) 



7r° yield extraction A 

Beam background A 

Tower energy scale B 

Tower gain spread B 

SMD energy scale B 

SMD gain spread B 

77/7r" B 

7r° yield fit B 



7.1 

1 (3) in d + Au 
3 
1 

12 
1 
2 

1.5 



in the BEMC could not be distinguished from genuine 
photons originating from the event vertex. Therefore, 
we varied the cutoff value for the electromagnetic energy 
fraction r in an event [Eq. ([T])] in the range r = 0.7-0.9. 
This propagated into 1-3% point-to-point systematic error 
of R^. 



VI. RESULTS AND DISCUSSION 
A. Cross section for neutral pion production 

The invariant differential cross section for tt" and rj 
production in inelastic p+ p interactions is given by 



E 



j3 P+P 
" "^incl 



dp3 



j3 P+P 



p+p 



2iTpT dpT dy 



(21) 



It has been shown that the singly diffractive contribu- 
tion to the inelastic cross section is negligible at pT > 
1 GeV/c [53]. Therefore, we can assume that the differen- 
tial inelastic cross section is equal to the differential NSD 
cross section in our px range. The total NSD cross section 
in p + p collisions was found to be c^sd = 30.0 ± 3.5 mb, 
and the total hadronic cross section in d -I- Au col lisions 
was found to be cr^'^^'^ = 2.21 ± 0.09 b (see section |lll a| . 

The measured cross sections for tt" production in the 
p + p (presented in Ref. [H] and included here for com- 
pleteness) and d + Au collisions are shown in Fig. 22 The 
cross sections are compared to the NLO pQCD calcu- 
lations [Sg. The CTEQ6M parton densities and 
the KKP fragmentation functions were used in the 
p+p calculation. The d + An calculation used the nuclear 
parton distributions for gold |65H67| . in addition. The 
factorization scale fi was set equal to pt and was varied by 
a factor of two to estimate the scale uncertainty, indicated 
by the dashed curves in the lower panels of Fig. [22j These 
panels show the ratio of the measured cross sections to the 
corresponding QCD predictions. The error bars shown in 
the plot are the statistical and the shaded bands are the 
systematic uncertainties. The normalization uncertain- 
ties are indicated by shaded bands around unity on the 
right-hand side of each ratio plot. The measured tt" cross 
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FIG. 22: (a) Cross section for inclusive tt^' production in 
p + p and d + Au collisions at ^sjvjv = 200 GeV, divided 
by the corresponding NLO pQCD calculations [52] for (b) 
p + p and (c) d + Au collisions. The error bars are statistical 
and shaded bands are pT-correlated systematic uncertainties. 
Normalization uncertainties are indicated by shaded bands 
around unity in the lower panels. 



sections were not corrected for feed-down contributions 



from T] — Stt , r] 



"7r°, and Kg -)■ 7r°7r° decays. 



which are expected to be negligible. It is seen that the 
measured tt" cross sections in both p + p and d + Au colli- 
sions are well described by the NLO pQCD calculations 
in the fragmentation region p^ > 2 GeV/c. 

In Fig. 23 we compare the tt*' measurements in the 
p + p and d + Au data with the previous tt^ measure- 
ments by STAR [23l [24| and with the 7r° measurements 
by PHENIX Here, and in all following figures, the 

cited data are shown with their statistical and systematic 
uncertainties added in quadrature. The normalization 
uncertainties shown by the grey bands in the figure are 
largely correlated between the 7r° and the tt* data points 
and uncorrelated with the PHENIX normalization uncer- 
tainties of similar magnitude. It is seen that the neutral 
and charged pion spectra from STAR agree very well 
in both p + p and d + Au data, in spite of different de- 
tector subsystems and analysis techniques used in these 
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FIG. 23: (color online) Cross section for inclusive tt" produc- 
tion in (a) p + p and (b) d + Au collisions at ^/snn = 200 GeV, 
divided by NLO pQCD calculations [52] and compared to the 
STAR 7r± I23l[24] and PHENIX 7r° [6l[26] measurements. The 
error bars are statistical and shaded bands are pT-correlated 
systematic uncertainties. Normalization uncertainties are in- 
dicated by shaded bands around unity in each panel. 
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FIG. 24: The J^/tt" ratio measured in (a) p + p and (b) d + Au 
collisions at ^/sNN = 200 GeV, compared to the PHENIX 
measurements [27 and to the niT scaling predictions. The 
error bars are statistical and shaded bands are pr-correlated 
systematic uncertainties. 



TABLE VI: The values of the power law fit parameters from 
Eq. (131 for the measured tt", rj, and 7inci cross sections. 



Data 



A (mb GeV-2 c^) po (GeV/c) 



n x^/ndf 



Tv°, p + p 
n°, d + An 

V, P + P 
rj, d + Au 

Tincl, P+P 
7incl 



(1.69±0.65)xl0^ 0.723±0.066 8.61±0.14 65/10 
(4.02±1.35)xl0^ 1.46±0.16 9.93±0.32 53/10 

(7.0±5.0)xl0i 1.33±0.23 9.83±0.44 30/9 
(3.33±0.41)xl0^ 1.33(fixed) 9.83(fixed) 32/10 
(3.1±0.1)xl0'' 0.941±0.268 8.61±0.40 2/5 
d + Au (2.4±0.1)xl0^ 0.697±0.126 7.88±0.23 2/5 



measurements. The present results extend the reach of 
STAR pion measurements to pr = 17 GeV/c. Compari- 
son to the cross sections measured by PHENIX shows 
good agreement, within errors, in both collision systems. 
However, we note that our data indicate a possible excess 
over the PHENIX measurements at pT > 10 GeV/c in 
both cases. 

To parametrize the pt dependence, the measured tt" 
cross section, as well as the ry and 7inci cross sections pre- 
sented in the following sections, were fitted to the power 
law function [Eq. ( |13[)] , and the resulting parameter values 
are listed in Table |Vll Because of the large uncertainties, 
the Pq and n parameters for the t] cross section in d + Au 
data had to be fixed at the corresponding p + p values to 
achieve a stable fit. The quoted values of x^/ndf indicate 
that these fits provide only a general guidance on the 
shapes of the spectra and do not necessarily describe all 
features seen in the data. 



In addition, the pure power law fit ^ Pj,"^ to the 7r° 
spectra at > 5 GeV/c gives m = 7.5 ± 0.1 (x^/ndf ~ 
6/5) iov p+p and m = 7.9 ± 0.2 (xV^df = 12/5) for 
d + Au collisions. 



B. Eta-to-pion ratio 



The 77 measurement is presented in Fig. 24 as the ratio 
of 77 to tt'^ invariant yields (shown in Figs. 16 and 15 
respectively). This allows many systematic uncertainties 
to cancel (see Table III I. The error definitions in the plot 



are the same as described above for the differential cross 
sections. The present measurement agrees very well with 
previous PHENIX results (open symbols) [57]. The solid 
lines show the asymptotic ratio R°° = 0.5, consistent 
with the world rj/n^ measurements (see Ref. [27] and 
references therein). The fit to our data for pT > 4: GeV/c 
gives Rn/^ = 0.46 ± 0.05 {p + p) and R^/^ = 0.44 ± 0.08 
(d+Au). The dashed curves in Fig.p4]show the prediction 
based on rriT scaling [25l [59H6T| . It is seen that the data 
are consistent with such scaling behavior. 



C. Nuclear modification factor 

A convenient way to observe medium-induced modifica- 
tion of particle production is to compare a nucleus-nucleus 
collision [A + B) with an incoherent superposition of the 
corresponding number of individual nucleon-nucleon colli- 
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FIG. 25: (color online) Nuclear modification factor R^a for (a) 
7v° and (b) ??, compared to the STAR 7r± [23|23 and PHENIX 
7r° measurements [271 . The error bars are statistical and 
shaded bands are pT-correlated systematic uncertainties. Nor- 
malization uncertainties are indicated by shaded bands around 
unity in each panel. 



sions {N + N). The nuclear modification factor Rab is 
defined as the ratio of the particle yield in nucleus-nucleus 
collisions and the yield in nucleon-nucleon collisions scaled 
with the number of binary collisions Acon, 



Rab 



(PNAB/dprdy 
{Tab) d^aP+P/dpTdy' 



(22) 



Here (Tab) is the nuclear overlap function, which is re- 
lated to the number of inelastic N + N collisions in one 
A + B collision through 



(Tab) cThZ = {Ncow) 



(23) 



In the absence of medium effects, the nuclear modification 
factor is unity, whereas Rab < 1 indicates a suppression 
of particle production in heavy-ion collisions, compared 
to an incoherent sum of nucleon-nucleon collisions. 
We calculated the RdA ratio [Eqs. (|22]) and (23l] as 



RdA 



<Zd^NaA/dpTdy 
{N,,n)d^aP+P/dpTdy' 



(24) 



where the nucleon-nucleon inelastic cross section was 
= 42 mb and {N^oii) = 7.5 ± 0.4 was 



taken to be CT;^^ 



calculated from the Glauber model (see section IIIC). 

The nuclear modification factors for tt" and 77 are shown 
in Fig. [25] The definition of the errors is the same as given 
for the differential cross sections in section [TVJI Figure [25] 
also shows the RdA for tt* measured by STAR PHI HI]- 
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FIG. 26: (color online) Nuclear modification factor Rcp mea- 
sured in d -I- Au collisions, compared to STAR tt^ measure- 
ment |24) . The error bars are statistical and shaded bands are 
PT-correlated systematic uncertainties. Common normaliza- 
tion uncertainty is indicated by a shaded band around unity. 



A good agreement between neutral and charged pion 
measurements by STAR is observed. Our tt*^ and 77 
data also agree reasonably well with the corresponding 
PHENIX measurements [?7l[55]. 

In peripheral d-t-Au collisions, the number of participant 
nucleons is small and the creation of a dense medium is not 
expected. This suggests that, instead otp+p interactions, 
peripheral collisions can be used as a reference. This was 
done through the ratio of particle production in 0-20% 
central (C) and 40-100% peripheral (P) events. 



R 



CP 



{N,,n)pd^Nc/dpTdy 
(^cou)c d^p/dprdy' 



(25) 



The advantage of this measure is that no p + p reference 
data are needed. The disadvantage is that a stronger 
model dependence is introduced due to the uncertainty in 
(A'coii)- Figure 26 shows the Rcp ratio for tt", compared 
to the STAR Tr^data [23]. It is seen that the agreement 
between the neutral and charged pion measurements is 
very good. The ratio stays constant at a value consistent 
with unity beyond pt = 8 GeV/c and, therefore, does not 
support a possible decrease of the ratio at high px, which 
was suggested by the tt* measurement. 



D. Direct photons 



The double ratio R^ [Eq. ( 15 )] measured in p + p and 
d + Au collisions is shown in Fig. [27] The shaded band 
near R^ = indicates our estimate of the upper limit of 
the remaining neutral hadron contamination. The curves 
correspond to NLO pQCD calculations jSHj, which were 
further evaluated as 



^'^Itheor ^ 



(7dir/7r°) 



NLO 



(7decay/7r°), 



(26) 



where the numerator is the ratio of the NLO pQCD direct 
photon and tt*' cross sections. The denominator is given 
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FIG. 27: The direct photon yield in (a) p + p and (b) d + Au 
collisions at = 200 GeV, expressed in terms of the 

double ratio i?^. The error bars are statistical and the shaded 
bands are pT-correlated systematic uncertainties. The curves 
correspond to NLO pQCD calculations of the differential 
cross sections for direct photon [B5] and tt" [55] production in 
p + p collisions for different factorization scales /j. (the upper 
pQCD curve corresponds to fi = pt/2). The upper limit of 
the fractional neutral hadron contamination Co is shown as 



the shaded band at R^, 



0. 



by the number of decay photons per tt^, as determined 



by the simulation described in section V E 



The NLO pQCD calculation used the CTEQ6M [53] 
parton densities and the GRV [51] parton-to-photon frag- 
mentation functions as an input. The scale dependence 
of this calculation, indicated by the dashed curves in the 
figure, was obtained by changing the scale fi in the calcula- 
tion of prompt photon production, while keeping the scale 
corresponding to the 7r° cross section fixed at /i — pt- In 
addition, we have varied the factorization scale for both 
cross sections simultaneously. The observed variation was 
quantitatively similar, although in the opposite direction. 
Since the measured 7r° spectrum favors the result of the 
pQCD calculation with /j, = pt, we have used this value 
for all three curves. 



Although Fig. 27 demonstrates that the measured val- 



ues of Rj are consistent with the calculated direct photon 
signal, the interpretation in this context has its limitations. 
First, the curves do not follow directly from the theory 
but depend on our simulation of the decay photon yields, 
as shown by Eq. (26 1. In addition, the NLO pQCD cross 
section for tt*^ production is less accurately constrained 
than that for prompt photon production. To allow for a 
more solid comparison to theoretical predictions, as well 
as to other experimental data, we have converted to 



an absolute cross section [Eq. (17)]. 

The calculation of absolute direct photon yields required 
that the systematic errors associated with inclusive photon 
yields, which canceled in the ratio R-y, were included again. 
We derived the 95% confidence limits for the cross section 



FIG. 28: The differential cross section for direct photon 
production at midrapidity in p -I- p and d + Au collisions 
at y/SNN = 200 GeV, compared to the PHENIX measure- 
ment [29] and to the NLO pQCD calculation |68j, which was 
scaled with (TdA) [Eq. ( |23[ )] in case of d + Au collisions. The 
error bars are statistical the shaded bands are pT-correlated 
systematic uncertainties. The arrows correspond to the 95% 
confidence limits, as defined in the text. 



in the pT bins where Rj did not correspond to a significant 
direct photon signal, assuming that the statistical and 
systematic errors both followed a Gaussian distribution, 
and using the fact that R^ > 1 by definition. 



Figure 28 shows the invariant cross section for direct 
photon production in p -\- p and d + Au collisions. The 
normalization uncertainties are not explicitly given in the 
figure. The NLO pQCD cross section for direct photon 
production in p + p collisions was scaled with the nuclear 
thickness function (T^a) [Eq. (23l] to account for the num- 
ber of binary collisions in the d-f- Au system. The precision 
of the presented measurement is limited by systematic 
uncertainties for pt < 9 GeV/c and by statistical uncer- 
tainties for larger px values. Nevertheless, our results are 
compatible with the NLO pQCD calculations. Our data 
are also in a good agreement with the direct photon cross 
section in p + p collisions measured by PHENIX [29] . 

Earlier measurements of direct photon production in 
proton-nucleus collisions have been performed by the 
E706 experiment [70j by scattering protons on a fixed 
beryllium target, with the proton beam energies of 530 
and 800 GeV. Those data show a strong discrepancy with 
pQCD calculations, which was attributed to multiple 
soft gluon radiation and phenomenologically described 
as an additional transverse impulse kx to the incoming 
partons [70] . It has also been argued that this discrepancy 
might be due to nuclear modifications present even in 
the light berillium nucleus used [7T] , although the p+p 
data by E706 show similar behavior at the same ^/snn ■ 
It is, therefore, of interest to compare our d + Au results 
to those of E706, as shown in Fig. [29] as a function of 
XT = 2pt/\/s, which is a suitable variable to compare 
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FIG. 30: The relative contribution of the fragmentation (Dj/g) 
and the pQCD hard scattering processes to the total direct 
photon cross section [68]. Both contributions are shown as a 
function of the pr of the produced photon. 



lead to a suppression of fragmentation photons, although 
an enhancement of directly produced photons has been 
proposed as well [731. Second, a measurement of the iden- 
tified prompt photons in p + p collisions is of interest for 
the RHIC spin program, a large part of which is devoted 
to constraining the gluon spin contribution to the spin 
of proton. The isolation criterion selects the quark-gluon 
Compton process and, therefore, enhances the sensitivity 
of the cross section to the gluon content of the proton. 



FIG. 29: (color online) (a) The cross sections per nucleon 
for direct photon production in d + An collisions, compared 
to the measurements in p + p and p + A collisions by E706 
experiment [70] at the comparable nucleon-nucleon center-of- 
mass energies. The theoretical curves were calculated with the 
INCNLO program ffjj with jj, = pt/2. All data have statistical 
and systematic uncertainties summed in quadrature. The ver- 
tical arrows indicate our estimate of the 95% confidence level, 
(b) The ratio of the data and the corresponding calculations. 



data taken with different beam energies. Whereas the 
ratio data/theory from Ref. [70] shows an increase of up 
to a factor of 4 towards low xt, our results at still lower 
Xt constrain such a potential deviation from theory to 
less than a factor of 2. It should be noted, however, that 
the data have been taken at significantly different ^/s. 

We have included both prompt and fragmentation com- 
ponents in the pQCD calculations, since our measurement 
was based on an inclusive sample of photons. The theo- 
retical calculation of these two components is shown in 



Fig. 30 A first measurement of the contribution from frag- 
mentation photons to the total direct photon cross section 
'm p + p collisions at RHIC was reported in Ref. 

The interest in disentangling photons from the frag- 
meirtation process and from the iiritial hard scattering 
is twofold. First, it has been observed that the hot and 
dense medium produced in central heavy-ion collisions at 
RHIC causes a suppression of particle yields, which has 
been attributed to induced gluon radiation from a par- 
ton traversing the medium. The same mechanism could 



E. Summary 

The present spectrum complements that of the tt^, 
which was measured by STAR in the transverse momen- 
tum range 0.35 < < 10 GcV/c, and extends up to 
Pt = 17 GeV/c. There is a good agreement between the 
neutral and charged pion cross sections in STAR, even 
though very different methods and detector subsystems 
were used. The 7r° cross section also agrees well with 
the measurements of PHENIX, and with the theoretical 
NLO pQCD calculations. 

This paper presents the first measurements of 77 me- 
son production by STAR, which are in agreement with 
the PHENIX measurements and with the niT scaling 
assumption. 

We present the measurements of the nuclear modifi- 
cation factor RdA, where the tt'^ production in d + Au 
collisions is compared to that m p + p, and Rcp, the com- 
parison between central and peripheral d + Au collisions. 
Both results are consistent with unity at high px, are in 
a good agreement with the tt* measurements, previously 
made by STAR, and significantly extend the pt range 
for light meson production measurements. 

This paper also reports the first measurement of direct 
photon production by STAR. A direct photon signal con- 
sistent with NLO pQCD calculation has been observed 
at high Pt for both systems. No strong modification of 
photon production in d -I- Au collisions was observed. 

The results will provide an important baseline for future 
Au + Au measurements in STAR. 



27 



Acknowledgments 

We thank the RHIC Operations Group and RCF at 

BNL, the NERSC Center at LBNL and the Open Sci- 
ence Grid consortium for providing resources and support. 
This work was supported in part by the Offices of NP 
and HEP within the U. S. DOE Office of Science, the U. S. 
NSF, the Sloan Foundation, the DFG cluster of excellence 
'Origin and Structure of the Universe', CNRS/IN2P3, 



STFC and EPSRC of the United Kingdom, FAPESP 
CNPq of Brazil, Ministry of Ed. and Sci. of the Rus- 
sian Federation, NNSFC, CAS, MoST, and MoE of 
China, GA and MSMT of the Czech Republic, FOM 
and NWO of the Netherlands, DAE, DST, and CSIR 
of India, Polish Ministry of Sci. and Higher Ed., Korea 
Research Foundation, Ministry of Sci., Ed. and Sports of 
the Rep. Of Croatia, Russian Ministry of Sci. and Tech, 
and RosAtom of Russia. 



[1] J. Adams et al. (STAR), Phys. Rev. Lett. 91, 172302 
(2003). 

[2] C. Adler et al. (STAR), Phys. Rev. Lett. 90, 082302 
(2003). 

[3] C. Adler et al. (STAR), Phys. Rev. Lett. 90, 032301 

(2003) . 

[4] J. Adams et al. (STAR), Phys. Rev. Lett. 93, 252301 

(2004) . 

[5] J. Adams et al. (STAR), Phys. Rev. Lett. 91, 072304 

(2003). 

[6] S. Adler et al. (PHENIX), Phys. Rev. Lett. 91, 072303 
(2003). 

[7] B. B. Back et al. (PHOBOS), Phys. Rev. Lett. 91, 072302 
(2003). 

[8] I. Arsene et al. (BRAHMS), Phys. Rev. Lett. 91, 072305 
(2003). 

[9] M. Gyulassy L Vitev, X.-N. Wang, and B.-W. Zhang, 
in Quark-Gluon Plasma, edited by R. C. Hwa and X.-N. 
Wang (2004), p. 123. 

[10] J. Cronin et al., Phys. Rev. D 11, 3105 (1975). 

[11] S. Y. Li and X.-N. Wang, Phys. Lett. B 527, 85 (2002). 

[12] K. Eskola, V. Kolhincn, and C. Salgado, Eur. Phys. J. 
C9, 61 (1999). 

[13] D. de Florian and R. Sassot, Phys. Rev. D 69 (2004). 
[14] B. Abelev et al. (STAR), Phys. Rev. C 76, 054903 (2007). 
[15] X.-N. Wang, Z. Huang, and I. Sarcevic, Phys. Rev. Lett. 
77, 231 (1996). 

[16] R. J. Fries, B. Miiller, and D. K. Srivastava, Phys. Rev. 

Lett. 90, 132301 (2003). 
[17] T. Peitzmann and M. H. Thoma, Phys. Rep. 364, 175 

(2002). 

[18] G. Bunce, N. Saito, J. Soffer, and W. Vogelsang, Ann. 

Rev. Nucl. Part. Sci. 50, 525 (2000). 
[19] B. Abelev et al. (STAR), Phys. Rev. D Rap. Comm. 80, 

111108 (2009). 

[20] J. Adams et al. (STAR), Phys. Rev. C 70, 044902 (2004). 
[21] B. I. Abelev et al. (STAR), Phys. Rev. C 80, 044905 
(2009). 

[22] J. Adams et al. (STAR), Phys. Rev. Lett. 97, 152302 

(2006) . 

[23] J. Adams et al. (STAR), Phys. Lett. B 616, 8 (2005). 
[24] J. Adams et al. (STAR), Phys. Lett. B 637, 161 (2006). 
[25] B. I. Abelev et al. (STAR), Phys. Rev. C 78, 044906 
(2008). 

[26] A. Adare et al. (PHENIX), Phys. Rev. D 76, 051106 

(2007) . 

[27] S. Adler et al. (PHENIX), Phys. Rev. C 75, 024909 (2007). 
[28] S. Adler et al. (PHENIX), Phys. Rev. Lett. 98, 172302 
(2007). 



[29] S. Adler et al. (PHENIX), Phys. Rev. Lett. 98, 012002 
(2007). 

[30] K. Ackermann et al. (STAR), Nucl. Instr. and Meth. 

A499, 624 (2003). 
[31] M. Anderson et al. (STAR), Nucl. Instr. and Meth. A499, 

659 (2003). 

[32] K. Ackermann et al. (STAR), Nucl. Instr. and Meth. 

A499, 713 (2003). 
[33] M. Beddo et al. (STAR), Nucl. Instr. and Meth. A499, 

725 (2003). 

[34] M. Bai et al., Nucl. Instr. and Meth. A499 (2003). 

[35] C. Adler et al., Nucl. Instr. and Meth. A470, 488 (2001). 

[36] J. Kiryluk (STAR), AIP Conf. Proc. 675, 424 (2003). 

[37] J. Kiryluk (STAR), 16th International Spin Physics Sym- 
posium Proc. pp. 718-721 (2005). 

[38] O. Grobenyuk, Ph.D. thesis, Utrecht University (2007), 
arXiv:0909.3006. 

[39] A. Drees and Z. Xu, Proceedings of the Particle Accelera- 
tor Conference p. 3120 (2001), Chicago, Illinois. 

[40] M. Miller, K. Reygers, S. Sanders, and P. Steinberg, Ann. 
Rev. Nucl. Part. Sci. 57, 205 (2007). 

[41] C. Adler et al. (STAR), Phys. Rev. Lett. 89, 202301 
(2002). 

[42] B. I. Abelev et al. (STAR), Phys. Rev. C 79, 034909 
(2009). 

[43] A. Baltz et al., Nucl. Instr. and Meth. A417, 1 (1998). 
[44] L. Hulthen and M. Sugawara, in Handhuch der Physik, 

vol. 39 (Springer- Verlag, Berlin, 1957). 
[45] R. Brun, R. Hagelberg, M. Hansroul, and J. Lassalle, 

Tech. Rep. CERN-DD-78-2-REV, CERN (1978). 
[46] J. Adams et al. (STAR), Phys. Rev. Lett. 92, 112301 

(2004). 

[47] T. Cormier, A. Pavlinov, M. Rykov, V. Rykov, and 
K. Shestermanov, Nucl. Instr. and Meth. A483, 734 

(2002). 

[48] A. Wetzler, Ph.D. thesis, J.W. Goethe-University (2006). 
[49] M. Russcher, Ph.D. thesis, Utrecht University (2008). 
[50] C. Amsler et al. (Particle Data Group), Phys. Lett. B 
667 (2008). 

[51] B. Abelev et al. (STAR), Phys. Rev. Lett. 97, 252001 

(2006). 

[52] B. .Jiigcr, A. Schafer, M. Stratmann, and W. Vogelsang, 

Phys. Rev. D 67, 054005 (2003). 
[53] T. Sjostrand, P. Eden, C. Priberg, L. Lonnblad, G. Miu, 

S. Mrenna, and E. Norrbin, Comput. Phys. Commun. 

135, 238 (2001). 
[54] G. D. Lafferty and T. R. Wyatt, Nucl. Instr. and Meth. 

A355, 541 (1995). 
[55] P. Aurenche et al., Phys. Rev. D 73, 094007 (2006). 



28 



[56] M. M. Aggarwal et al. (WA98), Phys. Rev. Lett. 85, 
003595 (2000). 

[57] S. Adler et al. (PHENIX), Phys. Rev. Lett. 94, 232301 
(2005). 

[58] A. Fasso, A. Ferrari, S. Roesler, P. R. Sala, F. Ballar- 
ini, A. Ottolenghi, G. Battistoni, F. Cerutti, E. Gadioli, 
M. V. Garzelli, et aL, Computing in High Energy and 
Nuclear Physics Conf. Proc. p. MOMT005 (2003), La 
JoUa, California. 

[59] M. Bourquin and J.-M. Gaillard, Nucl. Phys. B114, 334 
(1976). 

[60] J. Schaffner-Bielich, D. Kharzeev, L. McLerran, and 
R. Venugopalan, Nucl. Phys. A705, 494 (2002). 

[61] B. L Abelev et al. (STAR), Phys. Rev. C 75, 064901 
(2007). 

[62] S. Adler et al. (PHENIX), Phys. Rev. C 75, 051902 (2007). 
[63] J. Pumplin et al., J. High Energy Phys. 0207, 012 (2002). 
[64] B. Kniehl, G. Kramer, and B. Potter, Nucl. Phys. B852, 
514 (2000). 



[65] L. Frankfurt and M. Strikman, Eur. Phys. J. A5, 293 

(1999). 

[66] L. Frankfurt, V. Guzey, M. McDermott, and M. Strikman, 

J. High Energy Phys. 0202, 027 (2002). 
[67] L. Frankfurt, V. Guzey, and M. Strikman, Phys. Rev. D 

71, 054001 (2005). 
[68] L. E. Gordon and W. Vogelsang, Phys. Rev. D 48, 3136 

(1993). 

[69] M. Gliick, E. Rcya, and A. Vogt, Phys. Rev. D 48, 116 
(1993). 

[70] L. Apanasevich et al. (E706), Phys. Rev. D 70, 092009 
(2004). 

[71] G. Papp, P. Levai, and G. Fai, Phys. Rev. C 61, 021902 

(1999). 

[72] P. Aurenche, M. Fontannaz, J. P. Guillct, B. Kniehl, 
E. Pilon, and M. Werlen, Eur. Phys. J. C9, 107 (1999). 

[73] R. J. Fries and B. Miiller, Phys. Rev. Lett. 90, 132301 
(2003). 



